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Thermal  images  of  an  instrumented  ship  target  against  a  sea  background  have 
been  recorded  under  varying  environmental  conditions,  using  an  AGA  Thermovision 
780  radiometric  imaging  system  with  digital  data  recording.  Data  were  obtained  in  the 
8  -  14  Jim  spectral  band  at  an  angle  of  incidence  of  approximately  89  0  against  the  sea 
background.  These  pictures  have  been  analysed  together  with  measured  target 
emissivity  and  meteorological  parameters  to  give  the  temperature  and  radiance 
distributions  of  the  target.  The  influence  of  reflected  sky  radiation  on  the  background 
radiance  was  observed  under  varying  sea  surface  conditions.  A  predictive  model  of 
effective  radiance  temperature  difference  of  the  target  versus  sea  background  in  terms 
of  the  environmental  parameters  using  the  LOWTRAX  6  Radiance; Transmittance 
computer  code  was  developed.  This  effective  temperature  difference  has  been  shown  to 
vary  both  positively  and  negatively  compared  with  the  thermometric  temperature 
difference  depending  on  the  wind  speed  and  the  air  and  sea  surface  temperatures. 
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I.  INTRODUCTION 


Every  object  in  our  environment  emits  thermal  radiation.  The  peak  energy  of  this 
radiation  is  characteristic  of  the  object's  temperature  and,  for  temperatures  commonly 
found  in  our  environment,  is  concentrated  mostly  in  the  infrared  region  of  the 
electromagnetic  spectrum.  This  physical  phenomenon  has  found  many  applications 
since  the  second  world  war  as  it  enables  passive  detection  of  targets  by  then*  self- 
emission  alone. 

Passive  search,  detection,  identification  and  tracking  have  always  been  of  great 
importance  in  naval  warfare.  Also,  the  necessity  to  improve  our  combative  capabilities 
in  areas  where  radar  performance  is  degraded  (low  flying  small  cross  section  targets)  is 
imminent.  Infrared  technology  promises  much  in  this  area.  The  rapid  growth  of 
technology  in  solid  state  integrated  circuits,  microprocessors,  infrared  detectors  and  the 
better  understanding  of  infrared  physics  during  the  recent  decades  has  made  it  possible 
for  sophisticated  infrared  systems  to  be  designed  and  developed. 

A  very  important  step  in  designing  an  infrared  optical  system  is  the  knowledge  of 
the  spatial  and  spectral  radiances  of  the  target  and  its  background  as  well  as  the 
atmospheric  effects  on  the  emitted  thermal  radiation.  Almost  always,  targets  appear 
against  a  background  that  complicates  the  detection  process  since  radiation  from 
objects  surrounding  the  target  is  confused  with  radiation  emitted  or  reflected  from  the 
target.  In  the  analysis  of  thermal  imaging  systems  the  target-to-background  radiance 
contrast  is  frequently  described  by  the  equivalent  temperature  difference.  In  the  case  of 
observation  at  low  elevation  angles  above  the  sea  surface  the  surface  radiance  is 
strongly  influenced  by  the  reflection  of  the  radiation  emitted  by  the  sky.  The  actual 
thermometric  temperature  difference  is  no  longer  a  good  representation  of  the  radiance 
contrast.  An  "effective  temperature  difference"  can  however  be  defined  taking  into 
account  the  reflected  radiation.  Any  new  information  and  additional  data  about 
background  scenes  will  improve  the  understanding  of  these  effects  and  may  be  used  to 
optimize  computer  algorithms  or  to  implement  new  or  improved  techniques  in  signal 
processing  to  minimize  background  effects,  thus  improving  system  capabilities  and 
optimizing  performance. 


This  thesis  describes  the  collection  and  analysis  of  thermal  images  of  an 
instrumented  ship  target  and  its  background  under  various  well  defined  environmental 
conditions,  followed  by  the  development  of  a  model  for  the  effective  target-to- 
background  temperature  difference.  The  target  was  the  oceanographic  ship  R,  V 
"POINT  SUR"  and  the  background  was  sea  surface  and  sky.  The  data  were  collected, 
stored,  and  processed  using  the  "AGA  780  THERMOVISION",  a  liquid  nitrogen 
cooled  thermographic  device  with  digital  recording  and  data  analys;s  hardware.  All 
measurements  were  taken  from  a  distance  of  about  650  meters  using  the  Mercury 
Cadmium  Telluride  (HgCdTe)  detector  of  the  Thermovision  780  scanning  system  in  the 
8-14  pm  spectral  region.  Five  experiments  were  conducted  during  an  operational 
oceanography  student  cruise  on  mid  May  1986  and  an  adequate  number  of  thermal 
image  pictures  was  taken  for  various  em  ironmental  conditions  on  certain  dates. 

Analysing  these  thermal  images  using  the  measured  ship  surface  emissivitv  and 
the  AGA  software  package,  both  the  radiance  and  the  effective  temperature 
distributions  of  the  ship  were  obtained.  The  atmospheric  transmittance  for  this  purpose 
was  computed  using  the  LOWTRAN  6  code,  and  the  emisshity  was  measured  from 
samples  of  the  ship  paint.  The  effective  sea  surface  temperatuie  was  also  obtained  from 
these  pictures.  It  was  observed  that  the  apparent  sea  surface  temperature  was  affected 
by  reflected  sky  radiation,  and  a  model  was  developed  to  describe  the  effective  target- 
to-background  temperature  difference  under  varying  environmental  conditions. 

Later,  three  more  experiments  were  conducted  in  mid  November  19S6  in  order  to 
obtain  further,  mere  accurate,  data  under  different  environmental  conditions,  after  the 
Thermovision  had  been  calibrated.  This  was  done  to  check  the  previous  results  and 
improve  the  information  on  the  reflected  sky  thermal  radiation  and  the  effective  target- 
to-background  temperature  difference  relation. 

The  present  work  consists  of  eight  chapters.  In  the  first  chapter  we  state  the 
background  problem.  In  the  second  chapter  we  discuss  the  background  radiances  which 
affected  our  measurements,  the  marine  background  and  the  sky  background.  The 
atmospheric  propagation  problem  and  the  LOWTRAN  code  used  to  find  the 
atmospheric  transmittance  for  the  various  environmental  conditions  are  discussed  in 
chapter  three.  Chapter  four  provides  a  familiarization  with  the  equipment  used  to 
collect  the  data.  In  the  fifth  chapter  we  describe  the  calibration  problem  and  discuss 
the  derivation  of  the  empirical  calibration  relations.  In  the  sixth  chapter  we  discuss 
how  the  thermal  data  we  had  collected  were  analyzed  to  give  the  temperature  and 
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radiance  distribution  of  the  target.  A  further  analysis  of  these  d.ata,  to  see  how  the  sea 
background  radiance  was  affected  by  reflected  sky.  thermal  radiation,  is  given  in 
Chapter  :even.  In  this  chapter  we  also  calculate  the  apparent  >ea  surface  temperature 
based  on  environmental  parameters  and  compare  the  target-fo-background  temperature 
difference  sensed  by  our  detector  with  the  actual  thermometric  temperature  difference 
measurements.  In  the  last  chapter  rve  discuss  our  conclusions  and  recomendations  for 
further  work.  The  discussion  and  analysis  of  the  data  obtained  during  the  mid 
November  measurements  are  given  in  Appendix  A.  In  Appendix  B  the  TI-59  program 
used  to  get  the  in-band  radiant  flux  for  the  radiance  calculations  is  explained. 

The  results  of  this  work  will  lead  to  a  better  understanding  of  the  background 
effects  and  may  add  some  information  to  be  used  for  the  implementation  of  techniques 
in  signal  processing. 
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II.  TARGET  AND  BACKGROUND 


A.  BACKGROUND  RADIATION  NOISE 

Background  radiation  noise  is  noise  generated  in  an  optical  detector  by  photon 
flux  reaching  it  from  al1  points  other  than  file  target.  This  kind  of  noise  usually  is 
separated  into  two  components  [Ref.  Ij: 

1.,  Temporal,  which  is  due  to  fluctuations  in  the  background  emitted  photons. 

2.  Spatial,  which  is  due  to  background  structural  details. 

The  Temporal  background  noise  is  generated  by  the  fluctuation  in  the  number  of 
incident  photons  at  the  detector.  This  noise  usually  is  referred  to  as  background 
photon  noise  and  is  not  directly  related  to  the  target.  Detectors  sensitive  enough  to 
experience  background  photon  fluctuations  noise  are  called  BLIP  (  Background 
Limited  Photon  Detectors ). 

In  this  work  we  are  interested  in  the  spatial  component  of  the  background  noise, 
produced  by  the  structural  content  of  the  background  and  mainly  related  to  the  target. 
Therefore,  from  now  on  whenever  we  refer  to  background  noise  we  will  mean  the 
spatial  component. 

Since  a  target  is  always  surrounded  by  some  sort  of  background  it  is  obvious  that 
radiatipn  reaching  our  detector  from  this  background  is  confused  with  the  radiation 
emitted  or  reflected  from  the  target.  This  is  a  serious  problem  that  complicates  the 
detection  process  and  must  be  well  understood  for  designing  a  high  perforance  optical 
or  infrared  system.  [Ref.  1] 

Background  radiation  is  due  to  thermal  emission  and  reflected  or  scattered 
incident  radiation.  Many  parameters  affect  this  kind  of  radiation.  Of  high  importance 
among  them  are  the  changes  in  emissivity  and  reflectance  of  the  various  background 
surfaces  under  varied  environmental  conditions.  All  the  above  make  the  background 
problem  complicated  and  the  modeling  of  a  universal  background  function  impossible. 

Background  radiance  can  be  considered  mainly  as  two  categories: 

1.  Terrestrial  background 

2.  Celestial  background 

Terrestrial  background  refers  to  radiance  seen  by  infrared  systems  working  within 
the  atmosphere  of  the  earth.  This  kind  of  background  consists  of  radiation  that  several 
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surfaces  of  the  earth  are  emitting  or  reflecting  (  oceans,  clouds,  vegetation,  etc.  ).  It 
also  includes  sunlight  scattered  by  particles  in  the  atmosphere  and  atmospheric 
molecular  emission  (  sky  radiance  ). 

Celestial  refers  to  that  background  seen  by  a  system  working  in  space  and 
consists  of  radiation  from  the  stars  and  the  sun  reflectance  from  planets  and  space 
debris. 

Of  particular  interest  in  this  work  are  the  background  radiances  of  the  sea  surface 
and  the  sky.  Detailed  discussion  of  these  backgrounds  and  their  spectral  radiance 
characteristics,  especially  for  the  spectral  region  of  8  to  14  pm,  where  we  have  taken 
our  measurements,  is  given  in  Chapter  III. 

Since  target  and  background  radiances  travel  through  the  atmosphere  to  reach 
the  detector,  they  are  heavily  influenced  by  the  atmospheric  absorption  under  various 
meteorological  conditions.  Detailed  analysis  of  the  atmospheric  effects  on  thermal 
radiation  is  given  in  Chapter  IV. 


B.  BACKGROUND  RADIATION  CONTRAST 

The  target-to-background  temperature  difference  is  an  important  factor  in  the 
target  detection  probability  of  an  infrared  system.  The  higher  the  temperature 
difference  between  a  target  and  its  background,  the  greater 'the  distance  and  the 
probability  of  detection. 

This  temperature  difference  or  radiance  difference  is  usually  expressed  with  the 
quantity  Radiation  Contrast,  C^.  defined  in  [Ref  2]  as 

WT  -  WB 

CR  *  -  . 

wT  +  wB 

where  Wj,  and  Wg  are  the  target  (T)  and  background  (B)  radiant  emittances  (W/nr). 

The  8  to  14  pm  wavelength  radiation  contrast  versus  target-to-background 
temperature  difference  for  different  background  temperatures  is  shown  in  Fig  2.1  . 

When  the  temperature  difference  between  the  target  and  its  background  is  very 
small,  or  goes  to  zero,  the  radiation  contrast  goes  to  zero  and  the  detection  becomes 
very  difficult  if  not  impossible.  This  effect  is  known  as  the  washout  effect  and  happens 
for  "ambient  temperature"  targets  twice  a  day  due  to  the  solar  cycle. 
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Figure  2.1  Radiation  contrast  for  the  S-14  pm  region. 

In  [Ref.  3]  Hudson  describes  an  experiment  where  the  importance  of  the  contrast 
effect  is  shown.  In  that  experiment  a  vehicle  parked  in  an  open  field  was  viewed  by  an 
infrared  optical  system  over  a  period  of  24  hours.  During  that  observation,  contrast 
variations  were  found.  In  the  afternoon  the  contrast  was  positive,  since  the  vehicle 
heated  by  the  sun  was  warmer  than  its  background.  Later  in  the  early  evening  hours 
the  contrast  was  higher,  since  the  vehicle  with  larger  thermal  capacity  cooled  more 
slowly  than  the  background.  As  the  night  progressed,  the  contrast  passed  through  zero 
and  went  negative,  since  the  vehicle  continued  to  cool  rapidly  and  at  some  time  became 
colder  than  its  background.  During  the  morning  the  negative  contrast  was  still  present 
because  the  background  was  warming  more  rapidly.  After  some  hours  the  heating  of 
the  vehicle  was  enough  to  give  a  period  of  zero  contrast  and  then  turned  to  positive 
values  again. 

Hence,  twice  in  a  24  hour  period  we  are  lacking  sufficient  contrast  level  between 
the  target  and  its  background,  which  means  that  the  detection  of  any  target  at  these 
periods  is  not  possible,  and  since  this  effect  can  not  be  eliminated  we  have  to  utilize 
other  methods  of  detection. 

C.  BACKGROUND  SUPPRESSION 

Target  and  background  have  different  spectral  radiances  and  different  geometries, 
differences  that  make  it  easier  for  them  to  be  distinguished.  Whenever  a  spectral 
radiance  difference  exists  between  a  target  and  its  background  an  optical  spectral  filter 
can  be  used  to  reduce  the  unwanted  background  signals.  Even  if  spectral  filtering 
increases  the  target  to  background  ratio,  this  is  rarely  sufficient  to  render  a  system 
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operational  since  usually  the  irradiance  of  the  background  is  much  higher  than  that  of 
the  target.  A  common  characteristic  of  targets  is  that  they  are  always  smaller  in 
angular  extent  than  their  backgrounds,  a  typical  example  being  a  ship  against  the  sea. 
The  differences  in  geometry  between  targets  and  background  can  be  expressed  as 
differences  between  their  spatial  frequency  components.  For  such  geometry  differences, 
spatial  filtering  provides  a  good  suppression  of  backgrounds  of  this  extent  relative  to 
the  target  radiance. 

Spatial  filters  include  all  types  of  reticles  discussed  in  detail  in  [Ref.  3].  Spatial 
filtering  is  supplemented  by  spectral  filters  that  define  an  effective  spectral  bandpass 
for  the  system  and  electrical  filters  that  work  on  an  average  background  threshold  level 
to  chop  the  background  dc  signal. 

It  is  easily  understood  that  spatial  filtering  is  the  basic  suppression  scheme 
mainly  giving  discrimination  of  the  target  in  an  extended  background.  But  for  many 
military  systems,  depending  on  the  tasks  they  perform  (  identification,  tracking,  etc.), 
this  is  not  enough.  We  have  interference  from  other  smaller  kinds  of  background  such 
as  unwanted  objects  in  the  search  area,  for  example  buildings,  rocks,  clouds,  etc.  or 
unexpected  objects  crossing  our  detector's  field  of  view  such  as  birds,  different  kinds  of 
ships,  etc. 

Therefore,  more  sophisticated  and  complicated  spatial  filtering  systems  have  to 
be  used  for  all  the  different  purposes.  These  systems  are  supported  by  improved 
electronics  and  computer  software,  using  different  implementing  techniques  to  minimize 
the  background  effects  and  optimize  the  system's  performance. 
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hi.  Background  radiance 


a,  marine  background 

1.  Introduction 

The  sea  surface  infrared  radiance  (L$)  is  the  sum  of  the  radiance  caused  by  its 
thermal  emission  (Lw(l-R))  from  below  the  surface  and  the  radiance  due  to  reflected 
incident  radiation  (Lj.R)  (R:  reflectance  of  the  sea  surface). 


Ls  =  Lw(l-R)  +  L:.R 
Figure  3.1  Sea  surface  radiance. 

Major  factors  that  determine  the  marine  background  character  are  the 
following  [Ref.  4]: 

a.  Infrared  optical  properties  of  sea  water 

b.  Geometry  of  the  sea-surface  and  wave  slope  distribution 

c.  Sea  water  temperature  distribution 

d.  Properties  of  the  bottom  material . 

2.  Infrared  Optical  Properties  of  Sea  Water 

Thermal  radiation  from  the  sea  surface  occurs  at  wavelengths  having  their 
maxima  close  to  the  10  pm  region.  As  wc  can  see  in  Fig  3.2  taken  from  [Ref.  5].  at 
these  watelengths  water  is  essenfialiy  opaque  in  thickness  of  about  0.02  mm.  At 
was elengths  abo\e  3.5  fim  only  the  upper  lasers  of  the  sea,  of  thickness  about  0.0 1  cm 
and  less,  deternune  the  sea  surface  thermal  emission. 


REPRODUCED  AT  GOVERNMENT 


SAMPLE  THICKNESS  IN  MIC30NS 


Figure  j.2  Sea  water  aosorption  for  selected  wavelength  regions. 

The  radiation  emitted  from  the  sea  surface  is  dependent  on  the  underwater 
temperature  profile,  and  is  little  affected  by  surface  contamination.  Since  the  emissivitv 
is  equal  to  the  absorptivity  the  sea  surface  radiance  is  dependent  on  the  distribution  of 
the  absorption  coefficient  K  of  the  sea  water  as  shown  in  Fig  3.3  . 


Wavelength 


Figure  3.3  Absorption  coefficient  K  of  sea  water  versus  wavelength. 
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Subsurface  scattering  of  incident  radiation  is  not  significant  in  the  infrared 
range,  though  it  happens  in  the  solar  range  (0.35  -  3.0  pm).  In  that  region  the 
absorption  coefficient  has  a  very  low  value  and  radiation  penetrates  deeply  into  the 
water  and  is  scattered. 

Reflected  incident  radiation  is  composed  of  solar  radiation  and  sky  thermal 
radiation  .  Solar  radiation  dominates  at  short  wavelength  and  of  course  is  present  only 
in  the  daytime.  In  the  region  from  2  pm  to  4  pm  reflected  solar  emission  especially  at 
small  angles  of  incidence  is  greater  than  the  thermal  emission  of  the  sea  surface.  In 
comparison  to  this  region  only  about  0.0S  percent  of  the  sun  radiant  energy  is  found  in 
the  region  of  S  pm  to  14  pm  wavelengths,  which  means  that  only  reflected  sky  thermal 
radiation  influences  this  region  by  day  and  night. 

The  influence  of  reflected  sky  radiation  was  recognized  by  Me  Alister  [Ref.  6] 
who  found,  for  a  clear  sky  at  night,  an  apparent  sea  surface  temperature  a  few  tenths 
of  a  °C  colder  than  the  actual  temperature. 

Reflectance,  transmittance,  emissivity  and  index  of  refraction  for  water,  are 
shown  in  Fig  3.4  to  Fig  3.7  [Ref.  4].  Here  we  have  to  mention  that  distilled  water  and 
sea  water  transmittance  do  not  have  any  significant  difference. 

For  the  purpose  of  this  thesis,  where  we  have  worked  in  the  S  to  14  pm 
wavelength  region,  we  have  taken  into  account  only  reflected  sky  thermal  radiation. 


Figure  3.4  Transmittance  and  reflectivity  of  0.003  cm  of  sea  water. 
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REPRODUCED  AT  GOVERNMENT  EXPENSE 


Figure  3.5  Index  of  refraction  of  water. 


Figure  3.6  Reflection  from  water  surface  at  different  incidence  angles. 


Figure  3.7  Reflectance  and  emissivity  of  water. 
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3.  Geometry  of  the  Sea  Surface 

The  slope  distribution  of  the  waves  is  an  important  factor  in  the  reflection  and 
refraction  of  incident  radiation  on  a  roughened  sea  surface  at  various  wind  speeds.  In 
Fig  3.S,  taken  from  [Ref.  7],  is  shown  the  reflectance  of  solar  radiation  from  a  flat  sea 
surface  (  <y  =  0  )  and  from  seasurface  roughened  by  a  Beaufort  4  wind  (  a  =  0.2  ), 
where  cr  is  the  square  root  of  the  variance  of  the  wave  slope  distribution.  In  this  figure 
we  can  see  that  for  an  average  rough  sea  (  wind  speed  of  Beaufort  4  )  the  reflectance  of 
the  sea  near  the  horizon  goes  to  about  20%  where  as  the  reflectance  of  the  fiat  surface 
is  about  90%  and  above.  Consequently  we  can  see  how  high  the  emissivitv  (  1  -  R  )  of 
the  roughened  sea  surface  is  near  the  horizon.  Only  at  higher  angles  of  incidence,  about 
40°,  does  the  reflectance  of  the  roughened  sea  surface  go  slightly  higher  than  the 
reflectance  of  the  fiat  surface. 


Figure  3.8  Reflectance  from  a  fiat  (  cr  =  0  )  and  roughened  sea  surface  (cr  =  0.2  ). 
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In  [Ref.  7]  Charles  Cox  and  Walter  Munk  studied  this  phenomenon  for  solar 
radiation  in  the  visible  range.  They  assumed  unpolarized  sun  light  reflected  once  no 
double  scattering.  They  used  the  index  of  refraction  for  uncontaminated  sea  water  as  n 
=  1.338.  for  calculating  the  reflection  coefficient.  Information  is  lacking  for  the 
infrared  wavelength  region  and  for  night  conditions,  but  since  the  sky  radiance 
variation  is  almost  the  same  during  the  day  and  night,  and  the  reflectance  of  the  sea 
water  in  the  visible  region  of  electromagnetic  spectrum  at  zero  angle  of  incidence  is 
almost  equal  to  the  average  reflectance  in  the  infrared  region  from  2  pm  to  15  pm  (  see 
Fig  3.4  )  these  curves  are  considered  helpful.  For  the  purpose  of  this  thesis,  we  have 
used  these  curves  shown  in  Fig  3.8  from  [Ref.  7]  after  making  a  correction  for  the  index 
of  refraction  in  the  infrared  range  (  8  to  14  pm  )  we  have  used  in  our  work.  A  detailed 
discussion  of  this  is  given  in  Chapter  VIII. 

The  problem  with  these  curves  is  t’ .at  for  the  dashed  portions  near  the 
horizon,  above  S0°  angle  of  incidence,  the  function  being  used  is  not  strictly  applicable 
because  of  shadowing  and  multiple  reflections.  Also  the  lower  and  upper  branches  of 
the  (  cr  =  0.2  )  curve  represent  two  assumptions  used  in  [Ref.  7]  for  the  multiple 
reflection  effect.  The  basic  function  used  works  for  light  reflected  only  once  as  we 
mentioned  before.  True  values  are  expected  according  to  C.Cox  and  W.Munk  to  lie 
between  the  indicated  limits. 

The  curves  in  Fig  3.9  show  the  sea  surface  radiance  N(p)  divided  by  the  sky 
radiance  N$(0)  at  the  zenith,  due  to  sun  light  scatte  ng  by  air  molecules,  as  a  function 
of  the  vertical  angle  p,  for  a  clear  sky  and  an  overcast  sky,  and  is  also  taken  from 
[Ref.  7]. 

4.  Sea  Water  Temperature  Distribution 

The  sea  surface  temperature  determines  the  emitted  thermal  radiation  of  the 
sea.  This  temperature  is  about  0°C  (mainly  about  4°C)  in  arctic  regiuns  and  increases 
reaching  about  29  °C  near  the  equator.  Currents  may  produce  anomalies  of  some 
degree  in  temperature.  However  the  thermal  emission  of  the  sea  surface  has  its 
maximum  peak  near  the  10  pm  wavelength.  As  we  have  said  before  in  the  infrared 
optical  properties  section,  the  temperature  of  the  upper  layers  of  the  sea  surface,  of 
thickness  about  0.01  cm,  characterizes  the  thermal  emir.ance. 

Under  evaporation  conditions  the  temperature  of  the  upper  layer  of  0.1  mm  of 
the  sea  surface  was  measured  [Ref.  4]  as  colder  than  the  water  temperature  a  few 
centimeters  below  the  surface  of  about  0.6°C.  In  Fig  3.10  [Ref.  4]  are  shown 
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Figure  3.9  Sea  surface  radiance  divided  by  sky  radiance  at  zenith. 

measurements  of  typical  conditions  showing  the  sharpest  gradient  in  the  upper  layer  of 
1  mm. 

5.  Properties  of  the  Bottom  Materials 

*  • 

These  properties  are  important  in  regions  of  short  wavelengths  where 
radiation  can  penetrate  deep  into  the  water  causing  subsurface  scattering.  Such 
properties  are  not  important  for  the  work  done  in  this  thesis,  since  we  have  worked  in 
the  region  of  longer  wavelength  (  S  Jim  to  14  ftm  ). 
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Figure  3.10  Thermal  structure  of  sea  boundary  layer. 

B.  SPECTRAL  RADIANCE  OF  THE  SKY 
l.  Introduction 

Sky  background  radiance  in  the  infrared  is  caused  by  two  different  factors,  the 
scattering  of  radiation  from  the  sun  and  the  thermal  emission  of  the  .atmospheric 
constituents.  The  scattering  factor  is  important  in  the  short  wavelength  range,  close  to 
the  visible,  and  is  present  only  in  daytime.  The  thermal  emission  happens  at 
wavelengths  longer  than  4  pm  and  is  present  day  and  night. 

An  idealized  form  of  this  contribution  to  the  sky  spectral  radiance  is  given  in 
Fig  3.11  (Ref.  S].  This  figure  shows  the  spectral  radiance  of  a  biackbody  at  300°  K, 
close  to  the  temperature  of  the  earth,  and  the  spectral  radiance  of  a  biackbody  at 
6000°  K.  dose  to  the  temperature  of  the  sun.  On  estimating  the  sky  radiance  one  can 
assume  that  the  radiation  from  the  sun  is  diffusely  scattered  and  hence  the  sky  appears 
to  have  a  uniform  radiant  energy.  This  kind  of  maximum  brightness,  which  has  an 
average  radiance  about  2  xl0°  times  the  radiance  of  the  sun.  can  be  approximated  by 
a  bright  cloud  shown  as  the  lower  dashed  curve  in  the  figure.  But  normally  the  sky 
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does  not  scatter  this  large  amount  of  radiation;  a  nominal  value  is  about  one-  tenth  of 
this  maximum  spectral  radiance,  shown  as  the  solid  curve  on  -the  left,  named  "clear 
sky".  So  vve  can  describe  the  combined  effect  of  the  thermal  radiation  emitted  at  the 
temperature  of  the  earth  and  the  scattered  sunlight  as  combination  of  the  two  solid 
curves  shown  in  this  figure. 
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2.  Atmospheric  Emission 

a.  General 

The  atmospheric  emissivity  in  the  infrared  depends  on  the  amount  of  water 
vapor,  carbcn  dioxide  and  other  gases  which  give  rise  to  absorption  bands  due  to  their 
vibrationai  and  rotational  transitions.  Since  the  absorptivity  equals  the  emissivity 
(  Kirchhott's  hw  ).  we  can  predict  the  atmospheric  emission  from  the  absorption 
spectra  of  these  gases. 

b.  Spectral  Radiance  of  Char  Sky 

Sky  temperature  and  elevation  angle  are  the  main  factors  determining  the 
radiance  of  a  clear  sky.  The  elevation  angle  determines  the  length  of  the  atmospheiic 
path  and  hence  the  emissivity,  while  the  temperature  determines  the  blackbody 
emission  characteristics. 

The  effect  of  the  elevation  angle  on  spectral  radiance  is  presented  in  Fig 
3.12  [Ref.  $].  At  low  angles  of  elevation  (  near  the  horizon  )  we  can  see  nearly 
blackbody  radiation  at  the  ambient  temperature.  But  as  the  elevation  increases  above 
the  horizon  and  the  atmospheric  path  length  decreases  the  optical  density  decreases 
and  the  emissivity  is  no  longer  blackbody  (  c  =  1  )  everywhere.  In  the  regions 
centered  about  6.3  pm  (  water  vapor  absorption  band  )  and  about  15  pm  (  carbon 
dioxide  absorption  band  )  emissivity  is  very  high,  and  the  sky  remains  blackbody,  even 
at  the  zenith  (  90°).  But  between  these  regions,  in  the  region  from  about  8  pm  to  13 
pm  the  emissivity  is  poor  and  the  sky  is  no  longer  black.  The  general  background 
emission  m  this  region  is  due  to  the  wings  of  the  rvater  vapor  and  carbon  dioxide 
bands.  The  peak  that  we  see  at  about  9.6  pm  is  due  to  ozone  emission,  which 
increases  as  the  elevation  angle  decreases  from  the  zenith. 

Fig  3.12,  taken  from  [Ref.  8],  shows  spectra  of  clear  sky  radiance  measured 
in  September,  at  night,  from  a  height  of  11000  feet  above  sea  level,  at  the  ambient 
temperature  of  8°C.  The  angles  of  elevation  in  sequence  are  0°,  1.8°,  3.6°,  7.2°, 
14.2°,  30°  and  90°. 

The  effect  of  the  ambient  temperature  is  very  strong  for  this  kind  of 
radiation,  as  shown  in  Fig  3.13  .  This  figure  is  also  taken  from  [Ref.  8]  and  shows  the 
zenith  sky  spectral  radiance  for  two  different  ambient  temperatures.  The  dashed  curves 
are  the  spectral  blackbody  radiances  corresponding  to  the  ambient  temperatures  for  the 
two  separate  measurements  at  27.5  °C  and  2.5  °C.  These  measurements  were  taken  with 
different  amounts  of  air  mass,  which  affected  the  emissivity  and  hence  the  radiance  of 


Figure  3.12  Sky  radiance  for  several  elevation  angles  above  horizon. 

the  spectrum  in  the  regions  with  poor  emittance  [Ref.  S].  The  low  temperature 
measurement  was  taken  at  14110  feet  above  sea  level  and  the  higher  temperature 
measurement  at  6000  feet  above  sea  level.  The.  negative  values  shown  are  due  to 
measurement  error. 

c.  Spectral  Radiance  of  an  Overcast  Sky 

Overcast  skies  have  been  found  to  be  very  good  blackbody  radiators  (  z  = 

1  ),  having  a  temperature  at  low  altitudes  within  one  or  two  °C  from  the  ambient 
> 

surface  temperature  [Ref.  S]. 

As  mentioned  before,  water  in  the  infrared  is  a  good  absorber  (see  Fig  3.3) 
with  a  small  reflectivity.  Clouds  mostly  at  lower  altitudes  consist  of  water  droplets  with 
radii  from  1  pm  to  about  10  Jim  at  concentrations  above  300  cm0  with  temperatures 
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microns  wavelength 


Figure  3.13  Variation  of  zenith  sky  radiance  with  ambient  temperature. 

not  far  from  the  surface  temperature  [Ref.  9].  So  infrared  radiation  incident  on  a  cloud 
water  droplet  is  highly  absorbed  and  partially  reflected  and  transmitted.  However  due 
to  multiple  scattering  from  all  neighbouring  droplets  in  the  cloud  the  transmittance  of 
the  cloud  goes  essentially  to  zero. 

Hence  we  can  say  that  clouds  in  the  infrared  region  are  good  absorbers  and 
consequently  good  emitters,  causing  overcast  skies  to  behave  as  good  blackbody 
radiators. 
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d.  Spectral  Radiance  of  Clctuis 

Clouds  with  thickness  greater  than  100  m  and  clouds  with  water  droplet 
(  or  ice  crystal  for  higher  altitudes  )  counts  greater  than  0.5  g/cm~  do  not  transmit. 
Thus  absorptivity,  and  consequently  emissivity,  is  high,  approaching  one.  which  means 
that  clouds  show  a  very  good  blackbody  behavior  [Ref.  10]. 

Fig  3.14  [Ref.  10]  shows  the  emissivity  of  cloud  layers  versus  wavelength 
and  thickness  A Z  .  Also  in  Fig  3.15  [Ref.  10]  is  shown  how  the  integral  emissivity  (E), 
transmittance  (T)  and  absorptance  (A)  changes  with  thickness  AZ.  Here  we  can  clearly 
notice  that  at  thickness  AZ  above  60  m  the  transmittance  goes  to  zero  and  above  100 
m  the  emissivity  goes  to  zero. 

However,  looking  at  the  spectrum  of  the  under  side  of  a  cumulus  cloud  in  a 
clear  sky  (  Fig  3.16  )  at  an  elevation  angle  14.5°  [Ref.  3]  we  do  not  observe  a 
blackbody  spectrum  .  Since  the  emissivity  of  the  cloud  is  expected  to  be  high  and  since 
the  sky  is  partially  transparent  in  the  region  of  8  -  13pm  for  the  elevation  angle  of 
14.5°,  as  shown  before,  we  might  expect  to  see  clearly  the  emission  of  the  cloud  in  that 
spectral  region.  From  the  figure  we  see  that  the  emission  in  this  region  is  characteristic 
of  a  blackbody  at  -10°C  which  is  reasonable  for  the  temperature  of  the  cloud  (  lower 
dashed  line  in  figure  ).  What  is  really  happening  in  this  figure  is  that  we  have  a 
combination  of  the  atmospheric  emission  at  the  6  and  15  pm  wavelength  regions  from 
the  path  below  the  cloud  together  with  the  -10°  blackbody  emission  from  the  cloud  in 
the  wavelength  region  at  8  -  13  pm  . 

The  effects  of  thin  cirrus  clouds  in  a  clear  sky  modified  situation  is 
presented  in  Fig  3.18.  Here  we  can  see  that  the  emission  in  region  8-13  pm  has  been 
raised  in  comparison  with  the  clear  sky  situation  (  see  Fig  3.13  )  ,  and  also  the  ozone 
emission  is  reduced  since  it  comes  from  a  region  higher  than  the  clouds  and  is 
attenuated  passing  through  them. 

Fig  3.17,  from  [Ref.  8],  shows  spectra  measured  in  June  at  sea  level  with 
an  ambient  temperature  about  2S°C.  The  elevation  angles  are  0°,  1.8°,  3.6°,  7.2°. 
14.5°,  30°  and  90°. 

3.  Scattering  of  the  Skv  Radiation 

The  position  of  the  sun  relative  to  the  area  of  the  observed  sky  and  the  length 
of  the  scattering  path  due  to  the  elevation  angle  affect  strongly  the  amount  of  scattered 
radiation.  We  have  a  large  amount  of  scattering  close  to  the  sun  direction  as  well  as  an 
increase  in  scattered  radiance  with  a  decreasing  elevation  angle.  Decreasing  the 
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Figure  3.14  Spectral  emissivity  .of  cloud  layers  of  thickness  AZ. 


Figure  3.15  Integral  emissivity  (E). absorptivity  (A).transmittance  (T). 

elevation  angle  we  have  a  longer  atmospheric  path  from  which  scattering  can  be 

produced.  But  very  close  to  the  horizon  we  have  a  reduction  in  the  scattered  radiance 

mainly  because  of  water  vapor  absorption  which  reduces  the  effective  observation  path. 

This  kind  of  radiation  is  not  important  for  the  work  done  in  this  thesis  since  it 

occurs  at  short  wavelengths. 

*  w 
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Figure  3. It 


MICROWATTS  CM*  STCRADIAN4  MICRON*4 


IV.  ATMOSPHERIC  TRANSMITTANCE 


A.  ATMOSPHERIC  PROPAGATION 

Thermal  radiant  flux  from  a  target  has  to  pass  through  the  atmosphere  of  the 
earth  before  it  reaches  an  infrared  sensor.  During  this  passage,  the  thermal  radiation  of 
tne  target  is  attenuated  because  of  selective  absorption  due  to  several  atmospheric 
gases  and  scattering  by  small  particles  suspended  in  the  atrr.osphere.such  as  aerosols, 
rain,  snow,  smoke,  fog  and  haze.  Both  absorption  and  scattering  depend  on  wavelength 
(X). 

In  general  we  call  atmospheric  extinction  the  process  by  which  the  radiant  flux 
from  a  target  is  attenuated  passing  through  the  atmosphere.This  can  be  expressed  as: 

tA(X)  =  , 

where  t^(X)  is  the  atmospheric  transmission,  <J(X)  is  the  extinction  coefficient  and  R  is 
the  range.  The  extinction  coefficient  is  the  sum  of  the  absorption  coefiV-ent  «(X)  and 
the  scattering  coefficient  y(X) 

<r(X)  =  a(X)  +  y(X) . 

Both  absorption  and  scattering  cefficients  are  sums  of  molecular  and  aerosol 
components  [Ref.  2).  In  the  infrared  region  the  absorption  process  is  of  greater 
importance  (particularly  in  the  8  -  14  pm  )  than  the  scattering  process  which  is 
important  in  the  visible  and  near  infrared  region, 

Water  vapor,  carbon  dioxide  and  ozone  are  the  most  significant  absorbing  gases, 
having  broad  absorbing  bands  centered  at  the  specified  wavelengths:  water  \apor  (  2.7. 
3,2,  6.3  pm  ).  carbon  dioxide  (  2.7,  4.3,  15  pm  ),  ozone  (  4.8,  3.2,  6.3  pm  ).  The  bands 
of  carbon  dioxide  at  2.7  and  15  pm  and  water  vapor  at  6.3  pm  limit  the  atmospheric 
transmission  in  the  infrared  to  two  atmospheric  windows  at  3.5  to  5  pm  and  at  8  to  14 
pm. 

Infrared  systems  are  designed  to  operate  through  these  windows,  and  it  is  very 
important  for  the  performance  of  our  system  to  have  a  prediction  of  the  atmospheric 
attenuation  for  the  wavelength  and  the  environmental  conditions  in  which  we  are 
operating.  This  prediction  process  is  a  complicated  problem  and  has  been  faced  by 
computer  modeling  codes. 
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Such  models  are  the  high  resolution  models  (  HITRAN)  and  the  low  resolution 
model  (  LOWTRAN).  The  high  resolution  model  codes  provide  direct  computation 
with  high  resolution  single  frequency  results.  They  are  used  to  predict  atmospheric 
attenuation  for  laser  radiation.  The  low  resolution  model  codes  are  used  for  a  broad 
band  prediction  where  the  thermal  systems  work. 

LOWTRAN  code  was  used  to  predict  atmospheric  transmittance  for  the  work 
done  in  this  thesis  and  is  considered  in  the  next  sections  of  this  chapter. 

B.  LOWTRAN  6  CODE 

1.  General  Description 

LOWTRAN  6,  the  latest  edition  of  the  LOWTRAN  simulation  codes,  is  a 
FORTRAN  computer  .  curve-fitting  program,  relating  empirical  precomputed 
attenuation  data  to  simple  meteorological  and  environmental  inserted  parameters.  It 
was  designed  to  calculate  average  atmospheric  transmittance  and  thermal  radiance 
over  a  specified  broad  band,  ranging  from  359  to  40000  cm'*  (  0.25  to  28.57  Jim  )  with 
a  resolution  of  20  cm'*  averaging  bandwidth  in  steps  of  5  cm'*. 

The  average  transmittance  is  a  function  of  the  total  extinction,  the  product  of 
the  transmittances  due  to  water  vapor  line  absorption,  water  vapor  continuum 
absorption,  uniformly  mixed  gases  line  absorption,  nitrogen  continuum  absorption, 
aerosol  absorption,  aerosol  scattering  and  molecular  scattering.  A  consequence  of  the 
averaging  over  a  finite  spectral  band  is  that  the  rapid  fluctuations  of  the  atmospheric 
transmittance  in  that  band  vanish. 

Generally  we  can  say  that  the  LOWTRAN  code  is  a  fast  but  not  very 
accurate  computer  simulation  code. 

2.  LOWTRAN  6  Input  Data 

The  LOWTRAN  6  code  is  activated  by  submitting  four  cards  of  sequential 
data.  Each  card  provides  data  to  be  selected  and  asks  for  data  values  to  be  inserted  as 
follows: 

The  first  card  provides  a  selection  of  six  atmospheric  models.  These  models 
are  the  tropical  atmosphere,  midlatitude  summer,  midlatitude  winter,  subarctic  summer, 
subarctic  winter  and  the  1962  U.S.  standard  atmosphere.  Also  in  this  card  we  have  to 
indicate  the  type  of  atmospheric  path  (  horizontal,  vertical  or  slant  ),  and  determine  the 
execution  mode  (  transmittance  -  radiance  ). 
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The  second  card  provides  a  selection  of  the  type  of  extinction  for  nine 
bountary-  layer  aerosol  models  and  gives  a  default  meteorological  range.  These  models 
are 'the  rural  (vis  =  23  km),  rural  (vis  =  5  km),  navy  maritime  (computes  visibility), 
maritime  (vis  =  23  km),  urban  (vis  =  5  km),  tropospheric  (vis  =  50  km),  advection 
fog  (vis  =  0.2  km),  radiation  fog  (vis  =  0.5  km)  and  a  user  defined  model  (vis  =  23 
km).  Also  here  we  have  to  specify  the  air  mass  character  for  the  navy  maritime  model 
and  some  meteorological  conditions. 

The  third  card  defines  the  geometrical  path  and  we  have  to  insert  some 
parameters  like  initial  and  final  altitude,  initial  zenith  angle  and  length  of  our  path. 

In  the  fourth  card  we  have  to  specify  the  spectral  range  over  which  we  want 
data  for  atmospheric  transmittance  or  radiance  and  the  frequency  increment  over 
which  we  want  our  data  to  be  averaged. 

A  detailed  description  of  the  functions  of  these  cards  and  the  required  input 
data  is  given  in  [Ref.  11]. 

C.  LOWTRAN  6  CODE  USED  IN  THIS  WORK 

For  the  work  done  in  this  thesis  the  LOWTRAN  6  code  was  used  to  find  the 
atmospheric  transmittance  for  the  different  meteorological  conditions  we  experienced  . 

For  our  experimental  conditions,  mostly  overcast  skies  and  ambient  temperatures 
about  10°C  to  16°C  we  selected  the  midlatitude  winter  (  45°  N,  January  )  atmospheric 
model  for  horizontal  path.  Also  for  the  purpose  of  this  work  the  navy  maritime 
bountary  layer  aerosol  model  was  selected.  And  thus  because  this  model  takes  into 
account  the  existence  of  three  components  in  determining  the  aerosol  population.  First 
is  the  continental  component  which  is  sensitive  to  the  amount  of  time  the  air  mass  has 
spent  over  the  sea;  Second  is  the  stationary  component  which  depends  on  the  wind 
history  and  represents  aerosols  formed  by  high  winds  and  whitecap  conditions,  and  the 
third  one  is  the  fresh  component  which  is  characterised  by  the  current  wind  speed  and 
consists  of  whitecap  produced  droplets.  For  our  geometrical  path  we  inserted  target 
altitude  10  m  and  horizontal  path  length  650  m.  Our  spectral  range  was  710  -  1250 
cm'*  (8-14  pm  )  with  frequency  increment  5  cm'*. 

The  results  of  these  calculations  are  discussed  in  Chapter  VII. 
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V.  FAMILIARIZATION  WITH  THE  EQUIPMENT 


A.  GENERAL 

As  mentioned  in  the  introduction,  all  our  data  for  this  work  were  collected, 
stored  and  processed  using  the  AGA  Thermovision  780  system,  supplemented  with  a 
digital  microcomputer. 

The  system  is  designed  for  both  thermal  imaging  and  image  analysis.  It 
combines  real  infrared  scanning,  with  thermal  image  measurement  capabilities  and  is 
supplemented  by  a  microcomputer  for  digital  image  recording,  storing  and  processing. 

The  particular  system  we  used  was  a  dual  scanning  system  operating  in  the  3  - 
5.6  Jim  and  8  -  14  pm  windows  with  two  different  detectors.  For  the  purpose  of  this 
experimental  work  only  the  8-14  pm  scanner  was  used. 

B.  AGA  THERMOVISION  780 

1.  Basic  Description 

The  AGA  Thermovision  780  system  consists  of  a  real  time  scanner,  a  black 
and  white  monitor  chassis  display  and  a  digital  microcomputer. 

The  real  time  infrared  scanner  is  the  unit  that  converts  the  thermal  radiant 
energy  emitted  from  a  target  to  electronic  video  signals.  These  signals  are  amplified 
and  transmitted  through  a  cable  to  the  display  monitor  chassis  where  with  further 
amplification  they  produce  a  thermal  image  on  a  black  and  white  screen.  The  amplified 
signal  can  also  be  fed  to  a  color  monitor  through  an  interconnecting  device  in  order  to 
have  a  real  time  thermal  image  presentation.  Finally  video  and  data  signals  are 
transmitted  through  an  interconnecting  link  to  a  microcomputer.  Here  the  thermal 
image  is  desplayed  in  color  and  can  be  recorded,  stored  on  a  data  disk  and  processed 
for  temperature  evaluation  using  color  codes. 

2.  The  Infrared  Scanner 

There  are  two  kinds  of  scanning  unit  used  in  this  system,  the  short  wave  (SW) 
scanner  that  operates  in  the  3  to  5.6  pm  spectral  band  and  the  long  wave  (LW)  scanner 
that  operates  in  the  8  to  14  pm  band.  Scanner  systems  are  designed  with  single  or  dual 
scanner  units  which  cover  both  bands  simultaneously,  Our  scanner  as  we  have  already 
mentioned  was  a  dual  type.  Such  a  scanner  system  is  shown  in  Fig  5.1. 

Each  scanner  unit  consists  of  the  following  parts : 
a.  Electro-optical  scanning  mechanism 
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b.  Infrared  detector 

c.  Liquid  nitrogen  dewar  for  cooling  the  detector 

d.  Control  electronics  and  preamplifier 

A  simplified  block  diagram  of  the  scanner  unit  is  shown  in  Fig  5.2. 


Figure  5. 1  Dual  scanner  system- 


Figure  5.2  Simplified  scanner  unit  block  diagram. 
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Electromagnetic  thermal  energy  emitted  from  a  target  is  focussed  by  infrared 
lenses  into  a  vertical  prism.  The  lenses  for  the  short  wavelength  band  (SW)  are  made  of 
silicon  and  those  of  the  long  wavelength  band  (LW)  are  made  of  germanium.  Both 
kinds  are  covered  with  proper  antireflecting  coatings.  The  vertical  prism  is  rotated  at 
ISO  rpm  and  its  optical  output  is  passed  through  a  horizontal  prism  which  rotates  at 
1S000  rpm. 

Rotating  slotted  discs  control  the  rotation  of  both  prisms  and  are 
electronically  connected  to  the  horizontal  and  vertical  triggering  circuits  to  provide 
horizontal  and  vertical  triggering  pulses  for  the  monitor.  The  horizontal  and  vertical 
prisms  are  synchronized  in  such  a  way  that  four  fields  each  having  100  horizontal 
scanning  lines  produce  one  interlaced  frame.  Of  these  lines  only  70  lines  per  field  or 
280  per  frame  are  used  as  active  imaging  lines. 

The  output  from  the  horizontal  prism  passes  through  a  set  of  relay  optics 
containing  a  selectable  aperture  turret  with  eight  aperture  discs  between  f'1.8  and  f  20 
(fi'1.8  disc  was  selected  for  our  case  ).  It  then  passes  through  a  filter  unit  (  no  filter  for 
our  case  )  and  finally  is  focussed  on  a  single  detector  located  on  the  wall  of  the  dewar 
chamber  The  detectors  are  a  photovoltaic  Indium-Antimonide  (InSb)  for  the  short 
wavelength  band  and  a  photoconductive  Mercury  -  Cadmium  -  Telluride  (HgCdTe)  for 
the  long  wavelength  band.  Both  detectors  are  cooled  with  liquid  Nitrogen  which 
maintains  the  chamber  at  a  temperature  of -196° C. 

The  signal  output  from  the  detector  is  preamplified  within  the  scanner  unit 
and  the  video  signal  produced  is  applied  through  a  cable  to  the  black  and  white 
monitor  chassis  [Ref.  12]. 

The  arrangement  of  the  electro-optical  components  of  a  scanner  unit  is  shown 
in  Fig  5.3  [Ref.  12]. 

3.  Black  and  White  Monitor  Chassis 

The  video  output  signal  from  the  scanner  unit  is  applied  to  the  input 
connector  of  the  black  and  white  chassis.  This  signal  after  a  further  amplification  and 
processing  within  the  monitor  is  applied  to  the  display  screen. 

Outputs  from  this  monitor  chassis  are  used  for  a  color  monitor  connection 
and  the  interconnecting  link  for  the  microcomputer  connection  with  the  system. 

A  Black  and  White  monitor  chassis  is  shown  in  Fig  5.4  [Ref.  12]. 
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Figure  5.4  Black  and  White  monitor  chassis. 
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4.  Technical  Data 


Technical  data  for  the  AGA  Thermovision  780  were  found  in  [Ref.  13]  and  are 
presented  in  Table  1. 

5.  Digital  Image  Processing  System 

This  system  supplements  the  AGA  Thermovision  7S0  system  and  consis*  of 
an  interconnecting  link  and  the  microcomputer  which  utilizes  the  systems  software 
DISCO  3.0 

a.  Interconnecting  Link 

The  interconnecting  link  is  a  digital  on-line  interface  for  direct  interfacing 
between  the  Thermovision  780  and  the  IF  800  model  20  microcomputer.  The  main 
purpose  is  to  digitize  the  image,  hold  the  image  data  during  the  time  required  for 
transfer  from  the  Thermovision  to  the  microcomputer,  and  control  the-  data 
transferring  speed. 

b.  IF  800  Microcomputer 

IF  S00  is  a  business  level  personal  computer  which  uses  a  Z-S0A  ,  4  MHr 
clock  microprocessor  for  its  CPU.  It  comprises  a  keyboard,  a  printer,  a  display  unit 
and  a  mini  floppy  disk  drive. 

6.  DISCO  3.0  Software  Description 

The  software  used  by  this  system  consists  of  three  main  programs  the  DISCO, 
the  IMAGE  PROCESSING  and  the  UTILITY  programs.  Various  subprograms  as 
shown  in  the  DISCO  3.0  block  diagram  Fig  5.5  [Ref.  14]  supplement  the  main 
programs. 

The  main  program  DISCO  was  the  one  mostly  used  for  this  work  to  record, 
store  and  process  our  data.  Using  this  program  infrared  images  were  transferred  from 
the  Thermovision  7S0  system  to  the  microcomputer.  The  subprogram  FI  Record 
Image  was  used  to  record  these  thermal  images  and  display  them  in  color  on  the 
computer  screen.  Also  this  subprogram  was  used  to  enter  or  alter  our  parameters 
(emis«ivity,  atmospheric  transmittance  etc.),  to  select  and  enter  the  scanner 
identification  characteristics  used  and  store  the  images  on  a  data  disk.  The  subprogram 
F2  Evaluate  Image  was  used  to  evaluate  the  temperature  measurement  of  the  image 
using  color  scaling  with  tables  of  3  to  36  different  color  combinations.  Each  color  was 
assigned  a  particular  temperature  range.  The  standard  table  of  the  system  uses  8  colors 
and  was  the  one  used  in  our  case.  This  subprogram  was  also  used  to  display 
histograms  with  a  frequency  distribution  of  the  temperature  values,  to  display  single 


42 


REPRODUCED  AT  GOVERNMENT  EXPENSE 


TABLE  1 

AGA  THERMOVISION  7S0  TECHNICAL  CHARACTERISTICS 


FREFORMANCE 

Spectral  Range  i  3  -  S.  6  jlm  and/or  8-14  pm 
Frame  Ra'ce  s  6.25  sec"* 

Field  Rate  s  25  aOc"* 

Interlace  i  4  i  1 

REPLACEABLE  FORE  OFTICS 

FOV  Azimuth  *  i  ’• 

FOV  Elevation  *  i  7« 

1F0V  Azimuth  :  1.1  mrtd 

IFOV  Elevation  :  1.1  mrad 

HE.lO  i  0.  12*c  at  22*c 

MDT  i  <  0.  l«c 

Dynamic  Range  !  -20  to  900*c 

OFT I CAL  DATA 

Effective  Aperture  Area  :  24  cm2 

Aperture  Diameter  i  5.5  cm 

Effective  Focal  Length  s  9. 9  cm 

f /number  i  1.8 

DETECTOR 

Type  SW  : 

LW  : 

Humber  of  Elements  ! 

Feak  Wavelength  SK  i 

LW  : 

COOLING  SYSTEM 

Lguid  nitrogen  !  -196«c 

*  tlote:  The  system  can  take  five  different  fore  JenseL  with 
FOV  3. 5X3. 5,  7X7,  12X12,  20X20.  and  40X40.  The  l»ns 
we  us»d  was  the  one  with  FOV  7X7. 


Fhotovoltaic  InSb 
Fhotoconductl vj  HgCdTe 
1 

5  |im 
10  pm 
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DISCO  3.0 

Block  Diagram 


Block  diagram  of  DISCO  3.0  software. 


isotherms,  magnify  the  image,  transfer  a  new  image  on  the  display  from  those  which 
had  been  stored  on  the  data  disk  and  print  data  such  as  temperature  profiles  using  the 
system  printer.  For  the  last  case  our  system  had  some  limitations  since  the  system's 
maximum  capabilities  are  available  only  on  a  color  printer  which  was  not  available. 

C.  THERMAL  MEASUREMENT  TECHNIQUES 

1.  Introduction 

As  we  know  the  radiant  flux  density  received  by  an  infrared  system  does  not 
have  a  linear  relationship  to  the  object  temperature.  It  is  always  mixed  with  reflected  or 
emitted  background  radiation  and  is  attenuated  by  the  atmosphere.  Therefore,  certain 
calibration  and  correction  procedures  have  to  be  applied  to  a  particular  measurement 
in  order  to  have  an  accurate  temperature  evaluation. 

The  AGA  Thermovision  780  system  measures  thermal  radiation  between 
spectral  limits  as  mentioned  before.  The  numerical  value  of  the  received  radiation  is 
called  the  THERMAL  VALUE  and  is  measured  in  Isotherm  Units  (IU)  which  is  an 
arbitrary  unit  system.  This  thermal  value  is  proportional  to  the  received  photon  flux 
density,  but  is  not  proportional  to  the  temperature  of  the  target.  Therefore,  we  need  a 
calibration  function  as  a  graphical  curve  or  a  computer  program  to  relate  the  thermal 
value  (  isotherm  units )  with  the  temperature  of  the  target. 

2.  Calibration  Curves 

Calibration  curves  for  all  aperture  sizes  used  by  the  system  in  the  long 
wavelength  (  LW  )  spectral  band,  in  which  we  worked,  are  shown  in  Figs  5.6  through 
5.8  [Ref.  12].  These  curves  have  been  obtained  assuming  emissivity  t  =  1  (blackbody 
radiator  )  and  neglecting  atmospheric  damping  ta  =  1.  For  situations  different  from 
the  above  it  is  obvious  that  we  have  to  apply  correction  techniques. 

These  calibration  curves  are  accurately  described  by  the  following  equation 
taken  from  (Ref.  12] 

A 

I  =  -  ,  (eqn  5.1) 

Ce(S,'T).l 

where  I  is  the  thermal  value  in  isotherm  units  (  IU  )  for  the  absolute  temperature  (T) 
in  °K.  A,B.C  are  calibration  constants  that  depend  on  the  aperture,  filter,  scanner  etc. 

This  equation  was  used  by  the  microcomputer  program  to  calculate  the 
temperature  of  the  target  for  a  known  thermal  value  after  the  calibration  constants 
A.B,C  for  the  LW  shown  in  Table  2  had  been  inserted. 
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Figure  5.6  AGA  LW  calibration  curves  from  -20»C  to  50°C. 
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Figure  5.7  AG  A  LW  calibration  curves  from  -20  °C  to  150°C. 


ISOTHERM  UNITS 


Figure  5.8  AGA  LW  calibration  curves  from  0°C  to  1000°C. 


TABLE  2 


AGA  THERMOVISION  780 
LW  CALIBRATION  CURVE  CONSTANTS 


INDIVIDUAL  CALIBRATION  OF  780  DUAL  BBAR 


SERIAL  NUMBERS 

SCANNER  :  LWDR  4011 
DETECTOR  :  G  1739 
FILTER  :  NOF 

LENS  :  7  3105 

CALIBRATION  CONDITIONS 

AMBIENT  TEMPERATURE  :  24  C 

RELATIVE  HUMIDITY  :  55% 

OBJECT  DISTANCE  :  1.0  m 

CALIBRATION  DATE  :  82  -  09  -  27 


CALIBRATION  CURVE  CONSTANTS 


APERTURE 

A 

B 

C 

1.8 

-  3581 

1506.49 

-  0.436 

2.5 

-  4060 

1569.30 

-  0.759 

3.6 

-  6514 

1629.70 

-  1.824 

5.1 

5420 

1610.70 

2.796 

7.2 

1123 

1606.54 

1.09S 

10.0 

584 

1604.29 

0.885 

14.0 

306 

1610.60 

0.773 

20.0 

195 

1650.51 

0.767 

3.  Basic  Thermal  Measurements 

There  are  two  basic  methods  we  can  use  to  measure  temperature  using  the 
thermovision  780.  The  direct,  and  the  relative  measurement  method. 

In  the  direct  method  of  measurement  the  system  built  in  capabilities  (clamping 
and  temperature  compensation)  peimit  measurement  of  the  temperature  of  an  object 
without  using  any  other  reference  source  of  radiation. 

In  the  relative  measurement  method  the  radiation  of  an  external  reference 
source  with  known  temperature  and  emissivitv  is  compared  with  the  radiation  received 
from  the  object  to  obtain  the  object  temperature. 

The  direct  method  is  easier  to  use  but  not  as  accurate  as  the  relative  method, 
which,  however,  has  the  disadvantage  that  it  needs  an  accurate  reference  source  to  be 
available.  In  our  case  the  direct  method  was  used  since  we  didn't  have  any  reference 
source  with  known  temperature  and  emissivitv  located  near  the  target. 

4.  Direct  Measurement  Procedure 

In  this  subsection  we  will  describe  briefly  the  procedure  for  the  direct 
measurement  method,  using  the  Thermovision  first  alone  and  then  supplemented  with 
a  microcomputer  as  was  done  to  collect  data  for  this  work. 

To  aid  in  understanding  this  measurement  procedure,  the  monitor  chassis 
front  panel  controls  and  indicators  are  presented  in  Fig  5.9  [Ref.  12].  Subsequently  a 
short  description  of  the  operation  and  function  of  the  controls  and  indicators  referred 
to  in  this  measurement  procedure  is  given. 

a.  THERMAL  RANGE  (item  9  ):  This  is  a  nine  position  switch  calibrated 
between  2  and  1000  isotherm  units.  This  switch  selects  the  thermal  span  of 
interest. 

b.  THERMAL  LEVEL  (item  10  ):  This  control  sets  the  thermal  level  of  the  image. 

c.  PICTURE  MODE  "NOR?vIAL"  (  item  12b  ):  It  selects  the  normal  image 
displav,  where  the  cooler  object  temperature  appears  as  black  and  the  hotter  as 
white. 

d.  ISOTHERM  LEVEL  1  (item  13  ):  This  control  turns  on  the  isotherm  function 
and  identifies  the  discrete  thermal  levels  on  the  target.  The  isotherm  function  is 
indicated  by  saturated  white  areas,  which  represent  areas  of  the  same 
temperature.  A  marker  on  the  vertical  isotherm  scale  indicates  the  relative 
position  of  the  selected  level. 

e.  ISOTHERM  SCALE  (item  16  ):  This  is  a  scale  calibrated  from  -0.5  to  +0.5.  A 
reading  of  the  isotherm  marker  multiplied  by  the  thermal  range  gives  the 
relative  thermal  value  in  isotherm  units. 

For  both  cases  either  with  the  thermovision  itself  or  supplemented  by  the 
microcomputer,  starting  with  picture  mode  in  NORMAL,  we  have  to  adjust  the 
THERMAL  RANGE  and  LEVEL  controls  to  obtain  a  good  thermal  picture.  The 
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Figure  5.9  Monitor  chassis  front  panel  controls  and  indicators. 

picture  was  shown  on  the  Black/White  monitor  and  on  the  color  monitor  used  The 
measurement  procedure  is  shown  schematically  in  Fig  5.10  [Ref.  12].  First  with  the 
THERMAL  RANGE  control  switch  we  select  the  thermal  span  of  interest  in  isotherm 
units  (  2  to  1000  ),  about  the  thermal  level  of  the  image.  In  our  case  thermal  ranges  of 
2.  5  and  10  isotherm  units  were  selected.  Then  with  the  THERMAL  LEVEL  control 
knob  we  set  the  thermal  level  of  the  image.  At  the  point  at  which  a  satisfactory 
thermal  picture  is  obtained  the  corresponding  thermal  value  "L"  of  the  image  in 
isotherm  units,  is  measured  (i.e.,  the  thermal  value  of  the  isotherm  scale  at  zero  point, 
Fig  5.10). 

Using  the  thermovision  alone,  the  control  ISOTHERM  LEVEL  1  has  to  be 
* 

adjusted  to  brighten  up  the  point  of  interest  on  the  target  in  view.  A  marker  on  the 
vertical  scale  indicates  the  relative  position  of  the  selected  level.  The  reading  of  the 
isotherm  marker  multiplied  by  the  thermal  range  gives  the  relative  thermal  value  "i"  in 
isotherm  units,  shown  in  Fig  5.10.  Adding  the  two  values  "L"  and  "i"  we  have  the 


simplified  conditions,  having  =  1  and  £  =  1  ,  the  calibration  curve  can  be  used,  as 
in  Fig  5.10,  to  translate  the  measured  thermal  value  I'  to  temperature. 

Using  the  microcomputer  as  was  done  to  collect  data  for  this  w'ork  there  was 
no  need  to  go  through  the  isotherm  procedure.  The  computer  does  not  use  isotherms 
but  discrete  thermal  values.  Therefore,  we  had  the  isotherm  marker  at  zero.  The 
subprogram  FI  Record  Image  was  used  to  record  not  only  the  thermal  image  picture 
but  aFo  the  corresponding  thermal  value  and  thermal  range.  These  values  were  used  by 
the  subprogram  F2  Evaluate  Image  to  evaluate  and  present  on  the  display  the 
temperature  distribution  of  the  recorded  image  in  color  scaling,  using  the  standard 
eight  color  combination  table  of  the  system. 
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D.  COMPLEX  THERMAL  MEASUREMENTS 

1.  Introduction 

As  mentioned  before,  the  AGA  Thermovision  780  system  considers  an  ideal 
situation  for  the  values  it  measures.  But  as  we  know  there  are  several  factors  that 
influence  these  measurements.  So  the  true  object  measurement  must  be  derived  by  a 
calculation  in  which  we  must  take  into  account  all  the  factors  that  influence  the 
received  thermal  radiation. 

2.  Variables  Affecting  the  Measurement 

The  most  important  factors  that  affect  the. measured  thermal  radiation  are  the 
following: 

a.  Emissivity  of  the  target 

b.  Background  surrounding  the  target 

c.  Target  opacity 

d.  Target  size 

e.  Atmospheric  effects 

Most  of  the  objects  in  our  environment  are  not  blackbodies  and  hence  their 
emissivities  are  less  than  one.  For  our  case  the  emissivity  of  the  R/V  "POINT  SUR" 
was  found  to  be  about  0.95  as  will  be  explained  later.  The  influence  of  the  background 
was  discussed  in  detail  in  Chapters  11  and  III  and  the  atmospheric  effects  were 
discussed  in  Chapter  IV. 

The  target  opacity  has  to  be  taken  into  account  because  some  targets  are 
more  or  less  transparent  to  infrared  radiation  (  glass,  plastics,  clouds,  etc.  ).  Radiation 
received  from  objects  like  these  includes  radiation  transmitted  from  the  background. 
Also  the  emissivity  of  the  target  is  affected  and  is  hard  to  estimate.  For  our  case  the 
target  was  considered  as  an  opaque  Lambertian  surface. 

The  size  of  the  target  is  another  factor  to  be  taken  into  account.  Since  the 
detector  subtends  a  specific  solid  angle,  if  the  target  is  not  big  enough  to  cover  this 
angle  the  detector  receives  radiation  from  the  target  background.  This  will  tend  to 
affect  the  temperature  difference  indicated  between  the  object  and  the  background.  In 
our  case  the  target  was  not  big  enough  to  cover  the  field  of  view  of  a  7«X7«  lens  at  a 
distance  650  m.  For  most  measurements  it  occupied  1/3  to  2/5  of  the  display  screen. 

3.  Exact  Measurement  Formula  Derivation 

The  thermal  radiation  that  reaches  the  detectoi  of  our  system  is  the  sum  of 
the  target-emitted  thermal  radiation,  the  radiation  reflected  by  the  target  and  the 
radiation  emitted  from  the  atmosphere  as  shown  in  Fig  5.11. 


53 


Figure  5.1 1  Radiation  conditions  in  complex  measurement. 

If  we  assume  an  opaque  object  at  temperature  (TQ)  and  emissivity  (  eQ  ),  this 
object  emits  photon  flux  c0SQ.  SQ  is  the  photon  flux  of  a  blackbody  at  the  same  (  T0) 
temperature.  Since  the  above  flux  is  attenuated  by  atmospheric  transmittance  (  za  )  the 
object  radiation  reaching  our  detector  is  £0?aS0. 

Also  assuming  that  all  surrounding  surfaces  are  of  the  same  ambient 
temperature  (  Ta  )  which  determines  that  the  emissivity  is  close  to  one  (  ca  *  1  )  we 
can  say  that  the  radiation  reflected  by  the  object  is  P0Sa.  Sa  is  the  photon  flux  of  the 
surroundings  emitting  as  a  blackbody  (  ca  =  .1  )  at  temperature  (  Ta  )’and  p0  is  the 
reflectance  of  the  object.  From  the  first  assumption  that  the  object  is  opaque  we  get 
P0  =  1  -  e0  anc*  s*nce  ab°ve  radiation  is  attenuated  by  the  atmosphere  the 
reflected  radiation  reaching  our  detector  is  xa(  1  -  cQ  )Sa. 

Finally  the  atmospheric  radiation  is  given  by  (  1  -  ra  )Satm  where  Salm  is 
the  photon  flux  due  to  atmospheric  emission  and  scattering.  Hence  we  can  write  that 
the  photon  flux  S  Q  reaching  the  detector  is 

f  o  ~  eoT aso  +  Ta^  *  co^a  +  ^  '  ’a^atm  (ecin  • 

This  photon  flux  relation  can  be  converted  into  a  thermal  units  relation  since, 
as  we  mentioned  before,  the  thermal  value  (  IU  )  measured  by  the  AGA  Thermovision 
7S0  is  proportional  to  the  photon  flux  reaching  the  detector.  Therefore  we  can  write 
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the  relation  I  -  C  X  S  where  I  is  the  thermal  value,  C  the  proportionality  constant 
and  S  the  photon  flux. 

Substituting  S  =  I/C  in  Equation  5.1  we  have  the  same  equation  in  thermal 

values 


I'o  =  Wo  +  T«(l  -  £o)Ic£  +  -  T«)  him  5*3>  • 

This  equation  can  be  solved  for  IQ  to  find  the  exact  temperature  through  the 
calibration  curves  or  can  be  solved  to  find  the  emissivity  sQ  if  we  know  the  exact 
temperature  of  the  object  and  consequently  I  from  the  calibration  curves.  Therefore 
from  Equation  5.3  we  have 


I 


o 


V  0 


1 

(  T”  ’  1  ^atm  ’ 


(eqn  5.4)  or 


*  o  '  TcJa  *  "  Ta^atm 

£0  -  -  (eqn  5.6)  . 

Ta^o  ‘ 

For  our  case  the  computer  program  calculated  for  us  the  exact  temperature 
using  Equation  5.3  after  we  had  inserted  certain  values  for  the  parameters  emissivity, 
atmospheric  transmittance,  ambient  temperature  and  temperature  of  the  atmosphere. 


VI.  EMPIRICAL  CALIBRATION  OF  AGA  THERMOVISION  780 


A.  CALIBRATION  PROBLEM 

Trying  to  process  the  thermal  images  we  had  already  recorded,  we  realized  that 
the  temperatures  the  computer  was  calculating  were  much  higher  than  expected  for  the 
particular  scene. 

It  was  obvious  that  the  Thermovision  780  was  measuring  thermal  values  (  IU  ) 
higher  than  the  values  the  calibration  curves  were  giving  for  certain  target 
temperatures. 

To  solve  the  problem  we  had  to  derive  new  calibration  curves  to  use,  for  the 
exact  temperature  evaluation  of  the  stored  thermal  images.  At  that  time  we  had 
difficulty  deriving  new  calibration  constants  to  be  inserted  in  the  computer  program. 
Therefore,  we  tried  to  find  a  relation  between  the  exact  temperature  of  a  blackbody 
and  the  temperature  the  computer  program  was  calculating  using  the  faulty  calibration 
curves  for  the  blackbody  thermal  values  measurements.  In  this  way  we  had  an 
empirical  relation  as  a  curve  function,  accurate  enough  to  correct  the  temperatures  the 
computer  program  was  calculating  using  the  faulty  calibration  curve  on  processing  the 
thermal  images  we  had  stored. 

B.  THE  CALIBRATION  PROCESS 

1.  General 

For  the  calibration  procedure  we  required  a  blackbody  working  at  about  the 
temperatures  we  had  when  the  images  were  recorded.  The  air  temperature  at  that 
period  was  9.7  -  11.3°C,  the  sea  surface  temperature  was  11.2  -  I2.8°C  and  the 
laboratory  ambient  temperature,  where  the  calibration  procedure  was  carried  out,  was 
always  18  -  19  °  C.  Therefore  we  needed  to  have  a  blackbody  source  to  work  f  elow  the 
ambient  laboratory  temperature  in  order  to  get  the  calibration  relation  we  wanted. 

2.  Blackbody  Source 

For  the  purpose  of  this  work  wc  built  our  own  olackbodv  source  to  use  in  oa 
calibration  procedure.  As  we  have  already  mentioned  we  wanted  this  source  vo  work 
at  about  the  expected  target  temperatures,  which  were  below  the  Lboratorv  ambient 
temperature. 
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The  core  of  this  blackbody  source  was  a  cylinder  2S  cm  long  and  7.5  cm  wide 
made  of  thick  (  0.5  cm  )  aluminum.  One  end  of  this  cylinder  was  cut  at  an  angle  about 
60«  from  horizontal  and  it  was  covered  with  a  plate  of  the  same  thickness  of  aluminum 
for  minimizing  reflected  radiation  escaping  out  of  the  cavity.  The  other  end  was  cut  at 
90  o  and  it  was  covered  with  a  same  thickness  aluminum  plate  which  had  a  hole  of  2 
cm  diameter.  The  interior  was  painted  with  mat  high  emissivitv  black  paint.  Using  a 
nichrome  Mire  wound  uniformly  around  the  cylinder  we  could  heat  our  source  and  by 
pouring  liquid  nitrogen  inside  the  cylinder  we  could  cool  it  below  ambient  temperature. 
The  whole  source  was  covered  with  a  kind  of  fiber  glass  insulating  material  and  four 
thermocouples  were  used  to  check  the  temperature  stability  of  the  core  cylinder  and  the 
edge  plates.  In  Fig  6.1  a  schematic  representation  of  this  blackbody  source  is  shown. 

This  source  was  well  tested  showing  very  good  blackbody  behavior,  stabilizing 
and  maintaining  the  temperature  we  needed  quite  fast. 


Figure  6.1  Blcckbod}  source  used  for  calibration  procedure 
3.  Calibration 

We  wished  to  linJ  a  relation  between  the  exact  temperature  of  an  object  and 
the  temperature  calculated  by  the  computer  program  using  the  wrong  calibration 
curves  for  the  isothermal  values  measured  by  the  Thermovision  7S0. 


For  this  reason  our  blackbody  was  set  to  certain  temperatures  about  the 
expected  target  temperature  at  the  time  of  the  experiment.  After  these  temperatures 
were  stabilized  a  measurement  with  the  Thermovision  780  was  taken  from  a  distance  of 
about  2  m.  The  thermal  value  was  transferred  to  the  computer  and  the  thermal  image 
was  recorded,  to  be  processed  for  temperature  evaluation.  For  the  distance  of  2  m  the 
atmospheric  attenuation  was  neglected  (  ra  =  1  ).  The  above  procedure  was  repeated 
for  all  three  thermal  ranges  (  2.  5,  10  )  used  to  collect  data. 

The  above  calibration  temperature  measurements  for  the  three  different 
thermal  ranges  are  shown  in  the  following  Tables  3.4,5-.  Looking  at  these 
measurements  we  can  notice  that  there  is  a  linear  relation  between  the  blackbody  exact 
temperature  and  the  temperature  the  computer  was  calculating.  Using  the  method  of  a 
least  squares  linear  regression  we  found  the  equation  of  the  best  fitting  straight  line. 
This  line  was  the  calibration  curve  function  we  wanted  to  correct  the  temperatures  of 
the  stored  thermal  images. 

Plotted  calibration  data  measurements  and  the  best  fitting  straight  line  for  ail 
three  thermal  ranges  are  shown  in  the  following  Figures  6.2  to  6.4. 

C.  CONCLUSIONS 

From  the  calibration  curves  given  in  the  previous  chapter  we  can  see  that  the 
portion  of  the  cur’es  in  the  temperature  region  used-  for  the  calibration  procedure  is  a 
straight  line.  The  linearity  of  that  region  was  the  confirmation  that  the  linear  relation 
we  were  getting  was  reasonable. 

The  corrections  we  were  getting  from  these  straight  lines  were  accurate  to  within 
an  error  of  about  5%.  The  thermal  values  measured  by  the  thermovision  were 
accurate  to  within  half  a  thermal  unit.  This  caused  an  error  of  1.5  %  in  the  calculated 
temperatures.  The  variance  of  the  best  fitting  staight  line  for  each  thermal  range  caused 
a  3  %  error  in  the  thermal  range  2,  3.5  %  in  the  thermal  range  5  and  3.7  %  in  the 
thermal  range  10.  Therefore,  the  total  correction  error  was  within  4.5  %  in  the  thermal 
range  2,  5  %  in  the  thermal  range  5  and  5.2  %  in  the  thermal  range  10.  Of  course  we 
have  to  keep  in  mind  that  these  calibation  functions  are  only  valid  for  the  temperature 
region  about  our  target  temperatures  (  9  to  1S°C  ). 
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TABLE  3 


THERMAL  RANGE  2 

CALIBRATION  TEMPERATURE  MEASUREMENTS 


TEMPERATURES  IN  DEGREES  C  TEMPERATURES  IN  DEGREES  C 


BLACKBODY  AGA 
EXACT  COMPUTED 

DIFFERENCE 

BLACKBODY  AGA 
EXACT  COMPUTED 

DIFFERENCE 

3.0 

16.1 

13.1 

15.5 

26.9 

11.4 

4.2 

16.5 

12.3 

16.0 

27.3 

11.3 

5.0 

17.2 

12.2 

16.5 

27.6 

11.1 

5.6 

17.8 

12.2 

17.0 

28.0 

li.o 

6.3 

18.2 

11.9 

17.5 

2S.3 

10.S 

7.0 

18.7 

11.7 

18.0 

28.9 

10.9 

S.O 

19.7 

11.7 

18.5 

29.4  ' 

10.9 

8.5 

20.1 

11.6 

19.0 

29.8 

10.8 

9.0 

20.5 

11.3 

19.5 

30.2 

10.7 

9.5 

21.2 

11.6 

20.0 

30.7 

10.7 

10.0 

21.8 

11.8 

20.5 

31.1 

10.6 

11.0 

22.5 

11.5 

21.0 

31.6 

10.6 

12.0 

23.4 

11.5 

21.5 

32.0 

10.5 

12.5 

23.8 

11.3 

22.0 

32.4 

10.4 

13.0 

24.5 

11.5 

22.5 

33.0 

10.5 

13.5 

25.1 

11.6 

•  23.0 

33.6 

10.6 

14.0 

25.7 

11.7 

23.5 

33.9 

10.4 

14.5 

26.1 

11.6 

MEAN  DIFFERENCE  :  11.03 
STANDARD  DEVIATION  :  1.724 
VARIANCE  :  2.887 

BEST  FIT  STRAIGHT  LINE  :  Y  =  (1.0238)  X  -  11.62 
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Figure  6.2  Thermal  Range  2  calibration  temperature  measurements. 


TABLE  4 


THERMAL  RANGE  5 

CALIBRATION  TEMPERATURE  MEASUREMENTS 


TEMPERATURES  IN  DEGREES  C  TEMPERATURES  IN  DEGREES  C 


BLACKBODY 

AGA 

DIFFERENCE 

BLACKBODY 

AGA 

DIFFERENCE 

EXACT  COMPUTED 

EXACT  COMPUTED 

0.0 

14.1 

14.1 

16.0 

26.7 

10.7 

1.0 

15.1 

14.1 

17.0 

27.6 

10.6 

2.0 

15.8 

13.8 

18.0 

28.2 

10.2 

3.0 

16.8 

13.8 

19.0 

29.6 

10.6 

4.0 

17.6 

13.6 

20.0 

30.5 

10.5 

5.0 

18.2 

13.2 

21.0 

31.0 

10.0 

6.0 

18.9 

12.9 

22.0 

32.3 

10.3 

7.0 

19.2 

12.2 

'  23.0 

32.6 

9.6 

8.0 

20.4 

12.4 

24.0 

34.2 

10.2 

9.0 

21.7 

12.7 

25.0 

34.8 

9.8 

10.0 

22.3 

12.3 

26.0 

36.5 

10.5 

11.0 

23.0 

12.0 

27.0 

36.9 

9.9 

12.0 

24.2 

12.2 

28.0 

37.6 

9.6 

13.0 

25.0 

12.0 

29.0 

38.5 

9.; 

14.0 

25.7 

11.7 

MEAN  DIFFERENCE  :  11.55 
STANDARD  DEVIATION  :  1.531 
VARIANCE  :  2.265 

BEST  FIT  STRAIGHT  LINE  :  Y  -  (1.1956)  X  -  16.646 
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Thermal  Range  5  calibration  temperature  measurements. 


TABLE  5 


THERMAL  RANGE  10 

CALIBRATION  TEMPERATURE  MEASUREMETNS 


TEMPERATURES  IN  DEGREES  C  TEMPERATURES  IN  DEGREES  C 


BLACKBODY  AGA 

DIFFERENCE 

BLACKBODY  AGA 

DIFFERENCE 

EXACT 

COMPUTED 

EXACT 

COMPUTED 

0.5 

14.2 

•13.7 

16.5 

27.2 

10.7 

1.5 

15.0 

13.5 

17.5 

27.5 

10.0 

2.5 

15.6 

13.1 

19.5 

29.5 

10.0 

3.5 

16.8 

13.3 

20.5 

30.6 

10.1 

4.5 

17.7 

13.2 

21.5 

31.2 

9.7 

5.5 

17.9 

12.4 

22.5 

32.4 

9.9 

6.5 

18.6 

12.1 

23.5 

32.3 

9.3 

7.5 

19.4 

11.9 

24.5 

34.3 

9.3 

8.5 

20.1 

11.6 

25.5 

35.0 

9.5 

9.5 

20.6 

11.1 

26.5 

35.6 

9.1 

10.5 

21.0 

10.5 

27.5 

36.9 

9.4 

11.5 

22.9 

11.4 

28.5 

37.5 

9.0 

12.5 

24.1 

11.6 

29.5 

38.6 

9.1 

13.5 

25.2 

11.7 

14.5 

25.8 

11.3 

MEAN  DIFFERENCE  :  11.00 
STANDARD  DEVIATION  :  1.489 
VARIANCE  :  2.141 

BEST  FIT  STRAIGHT  LINE  :  Y  =  (1.1807)  X-  15.656 
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Figure  6.  < 


VII.  SIGNATURE  MEASUREMENTS 


A.  GENERAL 

All  signature  measurements  were  taken  during  mid  May  1986,  on  the  following 
dates  and  time  periods. 

1.  May  14  09:00  -  09:25 

2.  May  18  11:57-13:07 

3.  May  21  00:05  -  00:13 

4.  May  21  17:18-17:37 

5.  May  22  00:48  -  01:03 

Thermal  imaging  data  were  collected  for  the  spectral  range  8  -  14  pm  using  the 
LW  detector  of  the  AGA  Thermovision  780.  The  thermovision  system  and  its 
supporting  computer  were  set  up  in  a  laboratory  building  at  the  eastern  side  of  the 
HOPKINS  MARINE  STATION  grounds,  at  a  window  facing  directly  over  the  shore 
line  toward  the  closest  buoy.  The  laboratory  is  about  8  -  10  meters  above  sea  level. 

The  target  was  the  oceanographic  ship  R/V  POINT  SUR  which  was  on  an 
operational  oceanography  student  cruise.  During  those  periods  in  which  it  was  made 
available  for  signature  measurements  the  ship  maneuvered  in  the  vicinity  of  the  off 
shore  (  635  meters  )  buoy,  displaying  various  aspects  as  requested.  The  distance 
between  the  the  ship  and  our  detector  was  about  650  m.  At  that  distance,  the  closest 
the  ship  could  reach,  using  a  7°X7<>  FOVlens  we  had  an  image  covering  almost  2/3  of 
the  screen.  Six  thermocouples  installed  about  the  superstructure  were  used  to  have  a 
calibration  of  the  actual  ship  temperature.  The  thermocouples  were  sampled  at 
10-seconds  intervals  and  the  output  voltages  were  recorded,  to  be  later  converted  into 
temperatures. 

The  environmental  conditions  for  the  time  period  in  which  the  measurements 
were  taken  are  shown  in  the  data  presentation  tables.  Generally  for  these  periods  the 
area  was  under  the  influence  of  a  high  pressure  system,  with  North  West  winds  and 
coastal  fog,  especially  during  the  morning  hours. 

Thermal  images  for  different  ship  aspects,  were  recorded  and  stored  during  time 
periods  in  which  the  ship  was  available  for  us  and  the  environmental  conditions 
allowed.  All  data  were  recorded  assuming  the  ideal  conditions  appropriate  to  the 
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Thermovision  system  measurements,  that  the  target  was  a  biackbody  source  and 
atmospheric  attenuation  was  neglected. 

Fig  7.1  is  a  map  of  the-  Monterey  Peninsula  showing  the  area  where-  the 
signature  measurements  were  taken. 


B.  ‘  PARAMETERS  FOR  TEMPERATURE  EVALUATION 

As  mentioned  in  Chapter  five  the  processing  of  the  recorded  and  stored  thermal 
images  was  done  using  the  DISCO  3.0  computer  program.  For  temperature  evaluation 
the  program  asks  for  certain  parameters  to  be.  provided.  These  parameters,  as  we  have 
discussed  before  in  the  derivation  of  the  exact  measurement  formula,  in  Chapter  five, 
are  the  atmospheric  (air)  temperature,  ‘he  ambient  temperature,  the  atmospheric 
transmittance  and  the  target  emissivity. 

1.  Atmospheric  Transmittance 

The  atmospheric  transmittance  was  calculated  using  the  LOWTRAN  6  code. 
For  our  experimental  conditions,  temperatures  about  10  to  15°C,  we  selected  the 
midlafitude  winter  (  45°,  January  )  atmospheric  model  for  a  horizontal  path.  The  Navy 
Maritime  was  selected  for  the  boundary  layer  aerosol  model,  discussed  in  detail  in 
chapter  four.  For  the  geometrical  path  we  inserted  target  altitude  10  m  and  horizontal 
path  length  650  m.  The  spectral  range  was  set  at  8  *  14  um  (  710  -  1250  cm'*  )  with 
frequency  increment  5  cm'*. 

The  meteorological  data  needed,  mainly  wind  speed  and  wind  direction,  for 
determining  the  air  mass  character,  continental  inlluence,  were  taken  from  10-minute 
logged  meteorological  data  recorded  from  instruments  stationed  on  the  ship.  Such  data 
weren't  recorded  for  the  last  two  observation  periods  on  May  21  (  17:1S  -  17:37  )  and 
May  22  (  00:48  -  01:03  ).  For  these  periods  meteorological  conditions  were 
approximated  using  observations  taken  from  the  meteorological  station  at  POINT 
PINOS,  the  closest  station  in  the  area  of  the  experiment  (  see  Fig  7.1  ).  Since  the 
direction  of  the  wind  was  from  North  West  the  continental  influence  as  we  can  see 
from  the  map  in  Fig  7.1  was  not  very  strong.  So  on  a  scale  from  0  (  open  ocean  )  to  10 
(  maximum  continental  influence  )  we  inserted  mass  character  scale  from  2  to  3  . 

In  Table  6  are  shown  the  meteorological  data  used,  the  air  mass  character 
scale  (  ICSTL  )  and  the  calculated  atmospheric  tranmittance  for  each  observation 
period. 

2.  Emissivity 

The  emissivity  of  the  ship  surface  was  found  by  evaluating  the  emissivity  of 
two  metallic  pieces  taken  from  different  parts  on  the  ship.  These  specimens  were 
painted  with  slightly  different  white  paint,  representing  the  paint  coatings  of  the  ship 
surface.  These  had  also  been  expor"  i  under  the  same  environmental  conditions  for  the 
same  time  period  as  the  whole  ship  •  ce. 


67 


TABLE  6 


INSERTED  METEOROLOGICAL  DATA  AND  ATMOSPHERIC 

TRANSMITTANCE 

DATE  WIND  SPEED  WIND  DIR.  ICSTL  ATMOSP.  TRANSM. 


MAY  14  (09:00-09:25) 

4.80  m/s 

290  -  300 

-i 

0.S677 

MAY  18  (11:57-13:07) 

6.26  m/s 

270  -  280 

3 

0.8539 

MAY  21  (00:05-00:13) 

5.20  m/s 

270  -  300 

3 

0.8640 

MAY  21  (17:18-17:37) 

10.20  m/s 

320 

2 

0.S17S 

MAY  22  (00:48-01:03) 

9.25  m/s 

290  -  300 

** 

0.S282 

Assuming  that  the  emissivity  was  constant  over  the  spectral  region  involved, 
we  used  Equation  5.5  derived  in  Chapter  five,  for  our  calculations 

1  o  ‘  *  Ta^atm 

eo  = -  . 

Vo  '  V 

Making  some  more  assumptions  for  the  laboratory  environment  in  which  we  were 
working,  this  equation  was  simplified  further.  First  we  nad  to  assume  that  the  ambient 
temperature  and  the  atmospheric  temperature  were  the  same  in  our  ease.  This  meant 
that  the  equivalent  source  thermal  value  for  the  surroundings  (  Ia  )  had  to  be  the 
same  as  the  atmospheric  path  (  Iatm  ;.  So  the  above  equation  was  then  reduced  to 

*  o  *  *atm 
Eo 

Vo  *  V 

Second,  for  the  distance  of  about  2  m  between  the  detector  and  specimens  being 
measured,  we  assumed  no  atmospheric  attenuation,  so  tQ  =  1,  and  the  emissivity 
equation  was  reduced  again  to 
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where  I'0  is  the  Thermovision  780  measured  thermal  value,  Ia  is  the  calibrated 
thermal  value  for  the  known  ambient  temperature  and  IQ  is  the  calibrated  thermal 
value  for  the  known  object  temperature. 

With  this  last  equation  we  could  calculate  the  emissivity  of  the  specimens, 
accurately  enough,  if  we  knew  exactly  the  ambient  temperature  and  the  temperature  of 
the  specimens.  The  procedure  was  started  by  heating  the  specimens  as  uniformly  as  we 
could,  using  a  blower,  to  a  temperature  giving  good  thermal  images  above  the  ambient 
temperature.  This  temperature  was  checked  using  a  thermocouple.  With  the  best 
thermal  image  on  the  screen,  the  temperature  of  the  specimen,  the  ambient  temperature 
and  the  Thermovision  thermal  value  measurement  were  recorded.  The  ambient 
temperature  was  measured  using  a  common  thermometer. 

Knowing  the  exact  temperatures,  the  calibration  curves  with  an  appropriate 
correction,  discussed  in  -hapter  six  ,  were  used  to  get  the  thermal  values  (  IQ  )  for  the 
specimen  and  (  Ia  )  for  the  surroundings.  Therefore,  we  had  what  we  wanted  for  our 
calculation. 

The  above  procedure  was  repeated  several  times  at  different  temperatures  for 
each  specimen.  The  results  were  averaged  to  the  following  values. 

a.  0.943  first  specimen 

b.  0.958  second  specimen 

These  values  are  very  close  to  the  values  we  found  in  different  books  for  the 
emissivity  of  all  kinds  of  paint,  which  was  0.95.  For  the  calculations  done  in  this  work 
we  considered  the  emissivity  of  the  ship  surface  to  be  0.95,  taking  the  average  of  the 
two  values  we  found. 

3.  Atmospheric  Air  Temperature 

The  atmospheric  air  temperature  measurements  were  taken  from  logged 
meteorological  data  recorded  on  the  ship.  The  air  mass  temperatures  measured  close  to 
our  station  had  strong  continental  influence,  so  that  they  were  higher  than  the 
temperatures  measured  at  the  ship  location  and  not  representative  for  our  situation. 

The  a  mospheric  air  temperatures  are  shown  for  each  date  in  the  data 
presentation  taoles. 

4.  Ambient  Temperature 

These  temperatures  were  measured  using  a  common  mercury  thermometer 
near  the  open  window,  where  our  scanner  was  operating.  Their  value  was  between  the 
higher  room  temperature  and  the  lower  outside  temperature. 


69 


The  ambient  temperature  for  each  date  is  shown  in  the  data  presentation 

tables. 

C.  TEMPERATURE  DATA  EVALUATION  AND  PRESENTATION 

Inserting  the  correct  parameters,  the  computer  subprogram  Evaluate  image  was 
used  to  evaluate  the  temperature  distribution  of  the  ship,  for  each  thermal  picture 
recalled  from  storage.  This  temperature  distribution  was  presented  on  the  computer 
screen  with  the  standard  8  color  scaling  combination. 

For  better  presentation  of  these  temperatures  we  divided  the  ship  image  for  each 
aspect  in  small  squares  and  we  tried  to  assign  an  average  temperature  to  each  square, 
A  cartesian  system  was  used  to  present  these  temperatures  for  the  starboard  side,  the 
port  side  ard  the  stem.  These  temperatures  were  averaged  for  each  aspect,  for  each 
date,  in  order  to  have  the  final  temperature  distribution  for  the  particular  aspect  at  the 
time  period  our  data  were  averaged. 

The  problem  was  that  the  aspect  shown  by  the  ship  was  not  always  the  same, 
because  she  was  drifting  and  had  to  keep  in  station  using  the  engines  So  pictures,  say 
of  the  starboard  side,  taken  under  several  time  periods  had  different  aspect  angles  with 
the  line  of  sight.  This  caused  reflection  problems,  with  therma.  radiation  from  hot 
surfaces  such  as  the  stack  reflecting  from  cooler  surfaces.  Therefore,  our  scanner 
measured  not  only  the  emitted  thermal  radiation  but  both  the  reflected  and  the  emitted 
thermal  radiation  from  that  surface.  So  the  temperature  measurement  we  were 
evaluating  .vas  higher  than  the  correct  temperature  corresponding  to  that  surface.  To 
minimize  this  problem  we  tried  to  average  the  temperatures  only  from  those  pictures 
showing  the  most  perpendicular  aspect  to  the  line  of  sight. 

However  the  temperatures  obtained  were  not  correct,  as  discussed  in  detail  in 
Chapter  six,  since  the  computer-calculated  temperatures  had  still  to  be  corrected  using 
the  empirical  relation  already  derived.  Hence  the  average  temperature  for  each  square 
had  to  be  corrected  using  the  best  fit  straight  line  equation  corresponding  to  the 
thermal  range  in  which  the  data  were  collected. 

The  corrected  average  temperature  distribution  fox  the  starboard  side,  the  port 
side  and  the  smrn,  available  for  each  date,  are  given  in  the  following  Tables  7  to  11. 
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TABLE  7 

R/V  "POINT  SUR"  TEMPERATURE  MEASU REMET 
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TABLE  9 

R/V  "POINT  SUR"  TEMPERATURE  MEASUREMEN  TS 
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TABLE  10 

R/V  POINT  SUR"  TEMPERATURE  MEASUREMENTS  " 
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R/V  "POINT  SUR"  TEMPERATURE  MEASUREMENTS 
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D.  TEMPERATURE  DaTA  ACCURACY 

As  mentioned  before,  six  thermocouples  were  installed  about  the  superstuctiire, 
at  the  following  locations: 

T1  :  On  the  after  starboard  side  of  the  deckhouse 

T2  :  On  the  forward  starboard  side  of  the  deckhouse 

T3  :  On  the  aft  bulkhead  of  the  bridge  (  away  from  the  stack  ) 

T4  :  On  the  aft  bulkhead  of  the  deckhouse 
T5  :  On  the  port  side  of  the  deckhouse 
T6  :  On  the  after  deckhouse  forward  of  the  stack 
These  thermocouples  were  sampled  at  10-seconds  intervals  and  the  output 
voltages  were  recorded  and  later  converted  into  temperatures.  Thus,  we  had  an 
indication  of  the  actual  temperature  at  six  different  locations  on  the  ship,  for  the  whole 
time  period  we  were  taking  our  measurements.  In  this  way  we  had  an  absolute 
reference  to  compare  our  data  with  to  estimate  the  accuracy  of  the  process  of 
calculating  the  temperature  of  the  target  from  measurements  of  thermal  radiation. 

First  we  had  to  lit  the  thermocouples  in  the  squares  covering  the  areas  where 
they  were  located.  Then  we  had  to  find  the  mean  calculated  temperature  corresponding 
to  these  locations.  For  example  for  a  thermocouple  located  at  the  intersection  of  two 
or  four  squares,  the  calculated  temperature  corresponding  to  this  thermocouple  was 
taken  as  the  mean  temperature  of  the  two  or  four  squares.  Finally  we  had  to  compare 
the  above  temperature  with  the  thermocouple  temperature  averaged  over  for  the  same 
period  as  the  calculated  temperature  for  each  square. 

In  the  following  Table  12  are  shown  the  thermocouple  locations  and  the 
comparison  of  the  actual  with  the  measured  and  calculated  temperatures  in  degrees  C 
for  each  date. 

As  we  can  see,  for  most  of  the  time  the  temperature  differences  were  small.  The 
big  negative  differences  in  the  stern  area  for  the  two  first  dates  (  temperature  of 
thermocouple  higher  than  AGA  thermovision  temperature  ),  are  believed  to  occur 
because  the  thermocouples  were  not  yet  stabilized  at  that  period.  That  is,  at  a  time 
when  the  ship  was  turning,  hot  air  from  the  stack  area  was  blown  onto  these 
thermocouples,  making  their  measurements  fluctuate.  So  for  measurements  at  this  time, 
just  after  the  turn,  the  thermocouples  did  not  have  the  time  to  be  stabilized. 
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TABLE  12 
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E.  RADIANCE  CALCULATION 

Knowing  the  temperature  distribution  of  the  ship,  we  could  calculate  the  radiance 
for  the  spectral  region  of  the  measurements. 

Starting  with  the  Planck  blackbody  radiation  law  we  could  w'rite 


M  =  J  W(X,T)dX  =  |  (• 


2^hc 


,2 


X5 


-)  (- 


1 


aCh/XKT 


-)  dX 


-  1 


If  we  let  C  j  -  2jthc^  and  C2 
follows 


he 


K 


we  could  write  the  above  equation  as 


M-f 


W 


dX 


cm-1 


Solving  this  integral  in  the  spectral  region  8  -  14  pm.  for  the  known  temperatures 
we  found  the  in-band  radiant  flux  density  M  in  W/cnP\  With  the  assumption  that  our 
surface  was  an  opaque  Lambertian  surface,  we  could  divide  the  radiant  flux  by  ft  and 
multiply  ’by  the  emissivity  (  £  )  to  obtain  the  required  spectr?.l  radiance, 

CM  W 

L  =  -  - •  . 

2 

ft  enrsr 

The  integration  was  done  numerically  using  a  program  for  the  T1  -  59  calculator, 
details  of  which  are  g;ven  iu  Appendix  B.  The  results  of  our  calculations  are  shown  in 
the  following  Tables  13  to  17. 
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TABLE  14 

R/V  "POINT  SUR"  RADIANCE  MEASUREMENTS 
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TABLE  15 
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TABLE  17 
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1.93961  3.9537  14  3.9396  3  9813 


Vm.  DATA  ANALYSIS 


A.  INTRODUCTION 

In  this  chapter  we  win  analyze  how  the  sea  background  radiance  was  affected  due 
to  reflected  sky  thermal  radiation.  We  will  discuss  how  the  reflectance  and  the 
emissivity  of  a  sea  surface  roughened  by  various  wind  speeds  affected  the  background 
radiance.  Also  we  will  calculate  the  apparent  temperature  of  the  sea  surface,  and  will 
find  the  effective  target  to  background  temperature  difference  sensed  by  the 
Thermovision  under  the  specific  environmental  conditions  experienced. 

B.  EFFECTIVE  SEA  SURFACE  TEMPERATURE 

1.  Sea  Surface  Radiance 

As  mentioned  in  Chapter  ill  the  infnred  radiance  of  the  sea  is  the  sum  of  the 
radiance  due  to  thermal  emission  and  the  radiance  due  to  reflected  incident  radiation. 

The  reflected  incident  radiation  in  the  spectral  region  of  8  -  14  fim  is  only 
thermal  sky  radiation  by  day  and  night.  The  contribution  of  this  radiation  to  the  total 
sea  radiance  depends  on  the  reflectance  of  the  sea  surface  (Rj. 

Because  the  water  is  essentially  opaque  in  the  above  spectral  region  the 
reflectance  also  co  ^rols  the  emissivity  of  the  sea  surface,  since  for  an  opaque  material 
we  can  write  e  -  1  -  R  . 

Therefore,  we  can  write  that  the  total  radiance  of  the  sea  surface  is 

Ls  =  (  1  -  R  )  Lw  +  R  Lsky  ,  (eqn  8.1) 

where  Ls  is  the  radiance  of  the  sea  surface,  Lw  is  the  radiance  due  to  sea  water  thermal 
emission  and  Lskv  is  the  radiance  of  the  sky.  To  write  this  equation  we  assumed  that 
the  emissivity  of  the  sky,  as  we  have  previously  discussed  in  detail,  is  one. 

2.  Reflectance  of  the  Sea  Surface 

An  important  factor  in  dete  mming  the  reflectance  of  a  surface  is  its  geometry. 
For  the  sea  surface  this  geometry  depends  on  the  slope  distribution  of  the  waves  at 
various  wind  speeds. 

In  our  discussion  about  the  geometry  of  the  sea  surface  we  have  shown,  in  Fig 
3.8,  the  difference  in  reflectance  between  a  fiat  sea  surface  and  a  surface  roughened  by 
a  Beaufort  4  wind  speed. 
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These  resuits,  as  Fig  3.8,  were  obtained  by  Cox  and  Munk,  working  in  the 
visible  range  with  solar  radiation.  Such  information  is  lacking  for  the  infrared  and  night 
conditions,  although  the  above  curves  seem  to  be  helpful  in  this  area  also.  We  make 
this  assumption  since  the  reflectance  of  the  sea  water  at  zero  angle  of  incidence  is 
almost  equal  to  the  average  reflectance  in  the  infrared  region  from  2-15  pm  (  see  Fig 

3.4  ),  and  the  sky  radiance  variation  is  almost  the  same  for  day  and  night  conditions. 

For  our  analysis  we  used  the  same  curves  corrected  for  the  average  index  of 
refraction  in  the  spectral  region  8-14  fim.  Actually  ,T'e  wanted  to  find  a  correction 
factor  for  the  reflectance.  Since  we  knew  the  index  of  refraction  of  water  for  both  the 
visible  and  8  -  14  pm  regions,  we  could  find  the  corresponding  reflectances  at  normal 

,  9  9 

incidence  using  the  equation  R(0)  =  (  n  -  1  n  +  1  )  .  Their  ratio  gave  us  the 
required  correction  factor.  This  factor  was  used  to  correct  the  given  curves,  assuming 
that  the  correction  we  found  for  normal  incidence  applied  for  all  incidence  angles. 

[Ref.  7]  gives  the  index  of  refraction  used  to  calculate  the  reflectance,  for 
uncontaminated  sea  water,  as  n  =  1.338  and  quotes  the  ■  orresponding  calculated 
value  of  Ry(0)  as  0.020.  Even  though  the  reflectance  equation  for  n  =  1.338  gives 
Rv(0)  =  0.02089,  we  will  use  for  our  calculations  the  valueof  0.020  used  to  draw  the 
curves.  The  average  index  of  refraction  for  the  region  8  -  14  pm  was  found  from  Fig 

3.5  to  be  n  =  1.303  with  this  index  of  refraction  the  reflectance  was  calculated  as 
R-r(0)  =  0.01731.  As  we  can  see  the  reflectance  for  the  region  8  -  14  pm  at  normal 
incidence  is  lower  than  the  reflectance  in  the  visible  region,  in  the  ratio  of  Rir(0),  Ry(0) 
=  0.865.  With  this  correction  the  curves  were  shifted  a  little  lower. 

The  following  Tables  19  to  23  at  the  end  of  this  chapter  show  the  reflectance 
(R)  values  and  the  emissivity  (  1  -  R  )  values  as  they  were  found  from  the  corrected 
curves,  for  a  flat  sea  surface  and  one  roughened  by  a  Beaufort  4  wind,  for  certain 
angles  of  incidence.  Here  we  reiterate  that  the  valves  obtained  for  incidence  angles  near 
the  horizon  are  less  accurate,  because  the  curves  are  not  strictly  applicable  in  this  area, 
as  discussed  in  Chapter  HI. 

3.  Apparent  Temperature  Calculation 

As  we  said  before  in  section  one,  the  sea  surface  radiance  can  be  written  as 
the  following  sum  (  Eqn  8.1  ) 

L$  -  (  1  -  R  )  Lw  +  R  E sky  ’ 
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that  is,  that  the  radiant  flux  of  the  sea  surface  is  the  sum  of  the  emitted  radiant  flux 
and  the  reflected  incident  radiant  flux.  For  the  whole  spectral  region  we  can  write 
JTS4  =  (  1.-  R  )ffTw4  +  RcTsky4  ,  (eqn  8.2) 

where  T$  is  the  apparent  sea  surface  temperature,  Tw  is  the  actual  sea  surface 
temperature  and  Tskv  is  the  temperature  of  the  atmospheric  air.  Dividing  both  sides  of 
the  equation  by  c?  we  get  the  more  simplified  form 

Ts4  =  (  1  -  R  )  Tw4  +  R  Tsky4  (eqn  S.3)  . 

With  this  equation  we  could  calculate  the  apparent  sea  surface  temperature, 
the  temperature  that  our  thermal  system  would  sense  if  it  were  operating  in  the  whole 
spectrum  region  and  the  atmospheric  attenuation  was  neglected. 

With  the  above  assumption,  we  calculated  with  moderate  precision,  for  the 
region  of  observation,  the  apparent  sea  surface  temperature  for  a  flat  sea  surface  and 
one  roughened  by  Beaufort  4  wind.  Calculations  were  done  for  incident  angles  at  89°, 
where  we  had  been  working,  85°,  80°,  75°  and  70°.  The  reflectance  (R>  for  the  two 
kinds  of  surfaces  and  the  incidence  angles  listed  was  found  using  the  curves  in  Fig  3.8, 
after  being  corrected  for  the  spectral  region  8  -  14pm.  The  actual  sea  and  air 
temperatures  were  taken  from  the  recorded  meteorological  data  on  the  ship. 

The  results  of  these  calculations  are  shown  in  the  Tables  19  to  23  at  the  end 
of  the  chapter.  Detailed  analysis  and  discussion  of  those  results  is  given  in  the  following 
section. 

4.  Calculated  Temperature  Data  Discussion  and  Presentation 

From  the  equation  we  used  to  calculate  the  apparent  sea  surface  temperature 
we  can  see  that  the  two  terms  at  the  right  side  are  inversely  related  to  the  reflectance  of 
the  sea  surface. 

Therefore,  when  the  reflectance  of  the  sea  surface  is  high,  the  term  with  the 
sky  temperature  (  T$kv  )  dominates,  and  conversely  when  the  reflectance  is  low  the 
term  with  the  actual  sea  surface  temperature  (  T,v  )  dominates.  So  whenever  the 
reflectance  was  high  the  apparent  sea  surface  temperature  we  had  calculated  was  closer 
to  the  sky  temperature.  Since  for  our  case  the  sky  temperature  was  always  lower  than 
the  actual  sea  surface  temperature,  the  apparent  sea  surface  temperatures  we  have 
calculated  were  lower  than  the  actual.  On  the  contrary,  whenever  the  reflectance  was 
low,  which  means  that  the  emissivity  of  the  sea  surface  was  high,  the  apparent  sea 
surface  temperatures  we  have  calculated  were  close  to  the  actual  sea  surface 
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temperatures,  but  still  some  tenths  of  a  degree  lower,  because  of  the  small  contribution 
from  the  second  term. 

Looking  at  the  reflectance  values  given  for  the  two  kinds  of  surfaces,  at 
various  incident  angles,  in  Tables  19  to  23  we  see  that,  for  large  incidence  angles  (i.e., 
grading  incidence)  the  difference  in  reflectance  of  the  two  surfaces  is  large.  But  as  the 
incidence  angle  decreases  the  difference  decreases  and  at  about  70°  goes  to  zero,  (for 
more  details  see  Fig  3.8).  For  tnis  reason  the  difference  between  the  apparent 
temperatures  we  have  calculated  is  large  at  high  incidence  angles  (  89°,  85°  )  but  goes 
to  a  few  tenths  at  70  °. 

The  calculated  temperatures  are  shown  in  the  T$  row  of  the  Tables  19  to  23. 
In  the  T^  row  the  difference  between  the  actual  sea  surface  water  temperature  and  the 
calculated  sea  surface  temperature  is  shown.  Tm  is  the  temperature  computed  by  the 
Thermovision  780  system,  divided  by  the  atmospheric  attenuation.  This  was  done 
because  we  wanted  to  compare  this  temperature  with  the  T„  temperatures.  But  since 
the  Ts  temperatures  do  not  include  any  atmospheric  attenuation,  it  was  necessary' 
either  to  include  the  atmospheric  attenuation  into  these  values  or  amplify  the 
Thermovision  measurement,  to  obtain  a  value  without  the  atmospheric  attenuation. 
We  followed  the  second  method  and  we  got  the  values  Tm  shown  in  the  Tables  19  to 
23. 

Comparing  these  Tm  values  with  the  T$  values  we  have  calculated  for  incident 
angle  at  S9°,  we  can  see  clearly  how  much  the  reflectance  of  the  sea  surface, 
roughened  at  various  wind  speeds,  was  affecting  the  temperature  sensed  by  the 
Thermovision. 

Finally  in  the  Tm(j  row  of  the  Tables  19  to  23.  the  difference  between  the 
actual  sea  surface  temperature  and  the  Tm  temperature  is  shown.  As  we  have  expected 
this  difference  is  big  at  low  wind  "peeds  when  the  reflectance  is  high  and  the  appatent 
temperature  is  closer  to  the  sky  temperature. 

Generally  we  can  say  that  the  apparent  sea  surface  tempeiatures  found  were 
lower  than  the  actual  sea  surface  temperatures,  since  the  sky  was  always  colder.  Their 
difference  was  related  to  the  incidence  angle  and  the  roughness  of  the  sea  surface,  with 
higher  values  at  incidence  angles  near  the  horizon  and  lower  wind  speeds. 
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C.  TARGET-TO-BACKGROUND  TEMPERATURE  DIFFERENCE 
1.  Effective  Temperature  Difference 

An  important  factor  for  the  target  detection  probability  of  an  infrared  system, 
is  the  target-to-background  temperature  difference.  The  higher  the  temperature 
difference  between  a  target  and  its  background  the  greater  the  probability  and  the 
distance  of  detection. 

But  as  we  said  before  the  effective  temperature  difference  (  AT  )  that  our 
thermal  system  would  sense,  under  various  environmental  conditions,  looking  at 
various  incidence  angles,  might  not  be  the  same  as  the  actual  temperature  difference. 

Therefore,  the  precalculated  detection  distance,  using  the  actual  temperatures, 
at  which  we  expect  our  system  to  detect  a  specific  target,  may  not  coinside  with  the 
actual  distance  at  which  the  target  is  detected  by  the  system. 

In  this  analysis  we  attempted  to  find  the  effective  AT,  sensed  by  the 
Thermovision,  under  the  various  environmental  conditions  worked  in,  and  compare  it 
with  the  actual  AT  .  The  effective  AT  was  found  by  taking  the  difference  between  the 
temperatures  the  system  calculated  for  the  hottest  spot  on  the  ship,  ie.  the  stack,  and 
the  sea  surface.  These  temperatures,  were  calculated  under  the  ideal  conditions  in 
which  a  thermal  imaging  system  operates,  considering  emissivity  one  (  e  =  1  )  and 
neglecting  atmospheric  attenuation  (  xa  -  1  ).  and  were  corrected  using  the  calibration 
relation  we  had  derived.  The  actual  AT  was  found  using  the  temperature  of  the  stack, 
computed  from  the  radiation  measurements  using  the  exact  parameters  (emissivity, 
atmospheric  attenuation,  etc.  ). 

These  temperatures,  after  being  corrected  for  calibration  purposes,  were  tilt- 
most  accurate  available  for  our  case,  since  we  did  not  have  any  thermocouples  on  the 
stack  to  give  a  direct  indication  of  the  actual  temperature.  The  temperature  of  the  sea 
surface  was  taken  from  the  meteorological  data  recorded  on  the  ship. 

From  this  analysis  we  found  that  the  effective  Af,  for  the  environmental 
conditions  of  the  measurements  and  incidence  angles  near  the  horizon,  was  always 
higher,  by  about  a  factor  of  1.6  to  2  .than  the  actual  AT.  This  was  an  indication  that 
variations  in  emissivity  and  reflectance  of  the  sea  surface  under  the  experienced 
environmental  conditions  had  affected  the  effective  AT,  sensed  bv  tiie  Thermovision. 

Our  results  are  shown  in  Tables  19  to  23  at  the  end  of  the  chapter.  A  more 
detailed  discussion  of  these  results  is  given  in  the  following  section. 
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2.  Temperature  Difference  Data  Discussion  and  Presentation 

The  AGA  temperatures  shown  in  the  Tables  19  to  23,  for  the  hottest  spot  on 
the  ship  and  the  sea  surface  are  temperatures  calculated  by  the  computer  program 
using  the  thermal  values  measured  with  the  Thermovision  780,  assuming  £  =  1  and  r 
=  1  and  corrected  for  calibration. 

Therefore,  these  temperatures  correspond  to  the  thermal  radiation  attenuated 
by  the  atmospheric  effects  and  the  emissivity  of  the  measured  surface.  They  are  the 
effective  temperatures  our  system  was  sensing  under  the  environmental  conditions 
experienced  on  each  date,  at  high  incidence  angles  near  the  horizon  (  about  89  J  ) 

The  values  we  present  in  the  Tables  19  to  23  for  both  the  effective  and  actual 
temperatures  are  the  average  values  from  the  data  we  obtained  for  each  observation 
date. 

Comparing  the  effective  with  the  actual  AT  values  we  got  for  each  observation 
period,  we  see  that  the  effective  AT  for  the  environmental  condition:  experienced  was 
always  higher  than  the  actual.  It  is  obvious  that  there  are  two  main  factors  causing  the 
higher  effective  AT,  the  emissivity  difference  between  the  two  surfaces  and  the  reflected 
incident  radiation.  Since  the  emissivity  difference  was  small,  the  dominant  factor  was 
the  reflected  incident  radiation.  Because  the  reflectance  of  the  ship  surface  was  too 
small  (£  -  0.95)  to  cause  such  a  difference,  it  was  the  reflectance  of  the  sea  surface  we 
had  to  check.  But  we  have  already  done  this  in  the  previous  sections  the  calculatior  of 
the  apparent  sea  surface  temperature.  Then  the  conclusion  was  reached  that  the 
apparent  sea  surface  temperature  was  lower  than  the  actual  sea  surface  temperature 
because  the  sky  was  always  coider  than  the  sea  for  the  Leal  conditions.  These  lower 
temperatures  were  related  t  the  incidence  angle  and  the  roughness  of  the  sea  surface, 
with  the  lower  values  notice.,  at  incident  angles  near  the  horizon  and  low  wind  speeds. 

Hence  the  Thermovision  was  sensing  lower  sea  surface  temperatures  due  to 
reflected  thermal  radiation  from  the  colder  sky.  This  was  causing  the  effective  AT  to  be 
higher  than  the  actual.  Therefore,  we  would  e\p>-ct  these  lower  sea  surface 
temperatures  and  high  AT$  to  be  related  to  the  wind  speed  conditions  which  control 
the  roughness  and  consequently  the  reflectance  cf  the  sea  surface  since  tne  incidence 
angle  was  effectively  constant  (near  the  horizon,  about  89°  ). 

From  the  tabulated  and  plotted  (Fig  8.1  at  the  end  of  the  chapter)  ratio  of  die 
effective  and  the  actual  AT  for  the  various  w:od  speed  conditions  w?  experienced,  a 
significant  trend  can  be  seen.  Even  with  the  assumptions  and  cumulative  errors 
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inherent  in  the  approach  it  can  be  seen  that  for  the  lower  wind  speed  of  4.S  m,  s  oh  the 
14r^  of  May  we  had  an  effective  AT  higher  than  the  actual  by  about  a  factor  of  2.01 
and  for  the  higher  wind  speed  10.2  m/s  on  the  21st  of  May  (  i7:18  -  17:37  )  we  had  an 
effective  AT  higher  only  by  about  a  factor  1.63. 

Thus  as  expected  for  the  lower  wind  speed,  higher  reflectance,  we  had  the 
higher  AT  and  for  the  higher  wind  speed,  lower  reflectance,  lower  AT.  It  was  obvious 
that  the  effective  AT  was  strongly  dependent  on  the  reflected  incident  sky  radiation. 
Also  within  the  spread  of  the  plotted  points  we  may  conclude  that  the  AT  ratios  are 
proportional  to  the  wind  speeds. 

D.  CALCULATION  EXAMPLE 
1.  introduction 

To  help  understanding  of  the  following  Tables  19  to  23,  we  give  here  a 
calculation  example  tc  show  how  the  data  presented  were  obtained.  In  this  example  we 
use  the  information  we  had  for  the  18^  of  May  experiment. 

DATE  -  TIME  :  18th  MAY  1  i:57  -  13-07 

AIR  TEMPERATURE  (Tskv)  :  9.7°C 

SEA  SURFACE  TEMPERATURE  (Tw) :  13.1  °C 

WIND  SPEED  -  DIRECTION  :  6.26  m/s(B-4)  280° 

OVERCAST,  HAZY.  VISIBILITY  3  N'M 

ATMOSPHERIC  TRANSMITTANCE  :  0.8539 

Average  temperatures  computed  from  Thermovision  values  for  the  hottest 
point  on  the  ship  (Hot  Spot)  and  the  sea  surface  temperature  (S.S.T.),  (assuming  c  = 
i  and  -  1)  are  as  follows: 

UNCALIBRnTED  VALUES  :  HotSpot  :  22.S0°C 

S.S.T  :  20.65  °C 

CALIBRATED  VALUES  :  HotSpot  :  II.72°C 

S.S.T  :  9.52°C  . 

Aft  data  for  this  date  were  collected  using  THERMAL  RANGE  2  so  the  calibrated 
values  were  found  using  the  straight  line  equation  Y  =  (1.0238)  X  •  11.62  . 
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The  relation  between  the  Beaufort  sea  state  scale  and  wind  spee.d  (m,s)  is 
shown  below  in  Table  18. 


TABLE  IS 

BEAUFORT  SCALE 

B  WIND  SPEED  m/s 
0  0.00-  0.46 

1  0.51  -  1.54 

2  2.06  -  3.0S 

3  3.60  -  5.14 

4  5.65  -  8.23 

5  8.74  -  10.S0 

6  11.31-13.89 


2.  Apparent  Sea  Surface  Temperature 

To  calculate  the  sea  surface  apparent  temperature  we  used  Equation  8.3 

V  -  ( *  -  r  )  V  +  R  W> . 

Using  the  values  in  Fig  3.8  for  incidence  angle  89°,  where  we  have  been 
working,  we  got  the  reflectance  values  for  the  fiat  sea  surface  and  that  roughened  by  a 
Beaufort  4  wind.  These  values,  corrected  for  the  spectral  region  8  -  14  pm  using  the 
correction  factor  R-r/Rv  =  0.865,  are  the  following: 

FLAT  SURFACE  R  =  0.804 
ROUGHENED  SURFACE  R  -  0.110  . 

Therefore,  from  Equation  8.3  we  get  for  the  flat  surface 

Ts4  -  (  1  -  0.S04  )(  13.1  )4  +  0.110  (  9.7  )4 
Ts  -  10.65' °C 

and  for  the  rough  sea  surface  at  Beaufort  4  wind  speed  we  get 

T$4  -  (  1  -  0.110  )(  13.1  )4  +  0.1 10  (  9.7  )4 
T$  =  12.84°C. 

The  T^  values  represent  the  difference  (  Tw  -  T$  )  between  the  actual  sea 
surface  temperature  and  the  apparent  temperature  we  calculated.  Hence,  we  have  for 
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the  flat  surface  13.1  -  10.65  =  2.45  °C  and  13.1  -  12.84  =  0.26°  C  for  the  roughened 
surface. 

The  Tm  value  is  the  sea  surface  temperature  computed  by  the  AGA  system, 
but  without  the  atmospheric  attenuation.  To  get  this  value  we  proceeded  as  follows 
The  uncalibrated  value  which  the  system  calculated  for  the  sea  surface,  20.65  °C.  is 
amplified  by  the  atmospheric  transniission  tu  =  0.8539;  hence,  we  have 

20,65  *  0.8539  =  24.40°C  , 

This  value  was  corrected  for  calibration  purposes  ming  the  straight  line  equation  for 
thermal  range  2  ,  to  give 

Tm  =  24.40  (  1.0238  )  -  11.62  -  12.83°C  . 

As  we  would  expect,  this  value  must  be  close  to  the  Ts  value  for  the  roughened  surface 
(  B-4  ),  since  the  wind  speed  on  this  date  was  6.26  m/s  (  B-4  ). 

The  Tmcj  values  are  the  differences  Tw  -  Tm,  and  indicate  how  much  colder 
the  effective  temperature  sensed  by  the  Thermovision  was  than  the  actual  sea  surface 
temperature.  Hence  we  have 

Tmd  =  13.10  -  12.83  -  0.27°C  . 

3.  Effective  Sea  Surface  Temperature 

As  we  have  said  before,  the  values  given  for  the  AGA  are  the  values  the 
system  calculated,  corrected  for  calibration  .  These  values  are  not  corrected  for  the 
atmospheric  attenuation  and  emissivity.  Therefore,  for  the  hot  spot  we  have  the 
uncalibrated  value  22.8°C  -corrected  to  11.72°C,  using 

AGA  HOT  SPOT  =  22.8  {  1.0238)  -  11.62  =  11.72°C  ,. 

Similarly,  for  the  value  of  the  AGA  sea  surface  temperature  we  have  the  uncalibrated 
value  20.65 °C  which  corrected  gives  the  value  9.52 °C. 

For  the  actual  conditions,  the  temperature  of  the  hottest  spot  on  the  ship, 
since  we  did  not  have  any  thermocouples  on  this  location  ,  was  taken  as  the  highest 
temperature  we  found  for  this  location  in  the  calculated  temperature  distribution  of  the 
ship.  This  value  for  the  181^  of  May  was  14.28 °C.  The  sea  surface  temperature  was 
taken  from  the  recorded  data  on  the  ship. 


The  AT  value  represents  the  hot  spot  to  sea  surface  temperature  difference. 
Hence,  we  have  effective  ATe  =  2.20  6C,  from 

ATe  =  11.72  -  9.52  »  2.20°C  . 

The  actual  ATa  -  1.1S°C,  from 

ATa  =  14.2S- 13.10  «  1.1S°C  . 

The  ratio  ATe/ATa  =  1.864  represents  the  fractional  increase  in  effective  temperature 
for  the  particular  environ  nental  conditions  encounered. 
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TABLE  19 

SEA  SURFACE  TEMPERATURE  AND  EFFECTIVE  AT  "A' 


date  /  tike: 

WIND  SPEED  /  DIRECTION  : 
AIR  TEMPERATURE  : 
AMBIENT  TEMPERATURE  : 
SEA  SUREACE  TEMPERATURE  : 
OVERCAST,  HAZY,  VISIBILIT 
ATMOSPHERIC  TRANSMITTANCE 


14  MAY  09:00  -  09:25 

4.8  m/sec  (B-3)  290»-  300« 
10.  5»C 

15  »C 
12. 6»C 
1  4  m 

0. 8579 


.ANGLE 

OF  INCIDENCE  AND  SURFACE  CONDITIONS 

] 

II  89“ 

85“ 

80* 

75“ 

70“  | 

FLAT 

ESI 

8-4 

03312! 

3-4 

j  FLAT 

B-4 

R  Jj  0.804 

0.  110 

0.  513 

0.  138 

0.  292  io.  166 

0.  180 

0.  151 

o 

CS 

H 

6 

0.  120 

1-R  j |o. 196 

0.  390 

0.  487 

0.  862 

0.  820 

0.  849 

) 0. 880 

0.  880 

Ts  |  1 11.01 

12.42 

jll.  66 

12.36 

j  12.  09 1 12.  32 

1 12.  29 

12.35 

j- 12.  40 

12.  40 

Td  j  j  1.  59 

0.  18 

0.  94 

0.  24 1 

0.  5l|  0.28 

0.31 

0.  25 

j  0.20 

0.20 

l 


Tin 


11.  51 


Tmd 


1.  09 


TEMPER.  | 

HOT  SPOT 

SEA  SURFACE  . 

AT 

RATIO 

!  AGA  | 

IS.  32 

8.  44 

7.  38 

2.  01 

ACTUAL  j 

16.27 

12.  60 

3.  67 

DAT 

-4 

s 

d 

Tm 

Tmd 


REFLECTANCE  OP  THE  SEA  SURFACE 
FEAT  SEA  SURFACE 

SEA  SURFACE  ROUGHENED  BY  3EAUFORT  4  WIND 

CALCULATED  APPARENT  SEA  SURFACE  TEMPERATURE 

DIFFERENCE  BETWEEN  ACTUAL  AND  CALCULATED  SEA  TEMPERATURES 

AGA  SEA  TEMPERATURE  MEASUREMENT  WITHOUT  ATMOSPHERIC  ATTENUATION 

DIFFERENCE  3ETWEEN  ACTUAL  AND  AGA  SEA  TEMPERATURES 


*  ALL  TEMPERATURES  ARE  IN  °C 
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TABLE  20 

SEA  SURFACE  TEMPERATURE  AND  EFFECTIVE  AT  "B" 


DATE  /  TIME  :  13  MAY  11: 57  -  13: 07 

WIND  SPEED  /  DIRECTION  :  6.25  m/sec  (3-4)  270«-280» 

AIR  TEMPERATURE  :  9.  7»C 

AMBIENT  TEMPERATURE  :  17 *C 

SEA  SURFACE  TEMPERATURE  :  13. 1»C 

OVERCAST,  HAZY,  VISIBILITY  3  NM 

ATMOSPHERIC  TRANSMITTANCE  :  0.8539 


ANGLE  OF  INCIDENCE  AND  SURFACE  CONDITIONS 

■1 

89* 

•  as* 

80* 

75* 

70* 

■ 

Iflat 

3-4  | 

(flat 

3-4  j 

FLAT 

3-4 

FLAT 

3-4 

FLAT 

3-4 

KB 

m 

jo. 513 

0.  138 

0. 292  jo. 166 j 

0.  180 

0.  151 

0. 120 lo. 120 

l-R 

0.  196 

0. 890 j 

0.487  jo.  362 

jO.  708 

0.  334 

0.  320 !o. 849 

0.  380 l  o. 880 

Ts 

10.  55 

12.84 

11.  72 

12-  77  j 

12.37 

12.  70 

12.  60!  12.  74 

j 12. 32  12.32 

Td 

2.  45  j  0.26  j  1.  33 1  0.33 

0.73i  0.40 

0.44!  0.35 

•  0.28 

0.28 

Tn 

12.83 

Tmd 

0.27 

jfl  R 

TEMPER. 

HOT  SPOT 

SEA  SURFACE 

AT 

RATIO 

AGA 

11.  72 

9.52 

2.20 

1.  86 

ACTUAL 

14.28 

13.  10 

1.  18 

R  :  REFLECTANCE  OF  THE  SEA  SURFACE 
FLAT:  FLAT  SEA  SURFACE 

3-4  :  SEA  SURFACE  ROUGHENED  3Y  3EAUFORT  4  WIND. 

Ts  :  CALCULATED  APPARENT  SEA  SURFACE  TEMPERATURE 
Td  :  DIFFERENCE  BETWEEN  ACTUAL  AMD  CALCULATED  SEA  TEMPERATURES 
Tn  :  AGA  SEA  TEMPERATURE  MEASUREMENT  WITHOUT  ATMOSPHERIC  ATTENUATION 
Tmd  :  DIFFERENCE  3ETWEEN  ACTUAL  AND  AGA  SEA  TEMPERATURES 


*  ALL  TEMPERATURES  ARE  IN  °C 
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TABLE  21 

SEA  SURFACE  TEMPERATURE  AND  EFFECTS  E  AT  "C 


DATE  /  TIME  : 

WIND  SPEED  /  DIRECTION  : 
AIR  TEMPERATURE  : 

AMBIENT  TEMPERATURE  : 
SEA  SURFACE  TEMPERATURE  : 
CLEAR,  VISIBILITY  9  NM 
ATMOSPHERIC  TRANSMITTANCE 


21  MAY  00: OS  -  00:  13 

5.2  m/sec  (3-4)  270°-  300° 

11. 3°C 

17°C 

11. 9°C 

:  0.  8640 


'  ANGLE  OF  INCIDENCE  AND  SURFACE  CONDITIONS 


! 

j  89° 

00 

in 

o 

00 

o 

o 

O 

in 

70° 

L 

1 

FLAT 

FLAT 

B-4 

FLAT 

B-4 

FLAT 

3-4 

FLAT 

3-4 

R 

0.  804 

o.  no! 

0.  513 

0.  138 

C.  292 

0.  166 

0.  180 

0.  151 

0.  120 10.  120 

1-R 

1 

|0. 196 

0.  890 

0.  487 

0.  862 

0.  708 

0.  834 

1 

0.820 

0.  849 

_ 

0.  880 

0.  880 

i 

Ts  | 

11.  42 

11.  84 

mm 

11.  73 

11.  81 

11.  80 

11.81 

I 

11.  83 

11.  83 

Td 

j  0.48 

0.  06 

j  0.30 

0.  08 

0.  17  j  0.09 

0.  10 

0.  09 

i  0.07 

0.  07 

Tm 

11.  65 

Trad 

0.25 

[1 

TEMPER. 

HOT  SPOT 

SEA  SURFACE 

AT 

RATIO 

AGA 

12.  88 

7.  74 

5.14 

1.  90 

ACTUAL 

14.  60  j  11.  90 

2.  70 

1 

! 

R  :  REFLECTANCE  OF  THE  SEA  SURFACE 
FLAT:  FLAT  SEA  SURFACE 

3-4  :  SEA  SURFACE  ROUGHENED  BY  3EAUFORT  4  WIND 

Ts  :  CALCULATED  APPARENT  SEA  SURFACE  TEMPERATURE 

Td  :  DIFFERENCE  BETWEEN  ACTUAL  AND  CALCULATED  SEA  TEMPERATURES 

Tot  :  AGA  SEA  TEMPERATURE  MEASUREMENT  WITHOUT  ATMOSPHERIC  ATTENUATION 

Trad  :  DIFFERENCE  3ETWEEN  ACTUAL  AND  AGA  SEA  TEMPERATURES 


*  ALL  TEMPERATURES  ARE  IN  °C 
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TABLE  22 

SEA  SURFACE  TEMPERATURE  AND  EFFECTIVE  AT  "D" 


DATE  /  TIME 

WIND  SPEED  /  DIRECTION  : 
AIR  TEMPERATURE  : 
AMBIENT  TEMPERATURE  : 
SEA  SURFACE  TEMPERATURE  : 
SUN,  CLEAR,  VISIBILITY  10 
ATMOSPHERIC  TRANSMITTANCE 


21  MAY  17: 18  -  17: 37 
10.2  m/sec  (B-S)  320° 

12.  2°C 
19°C 

13.  2°C 
NM 

:  0. 8178 


ANGLE 

OF  INCIDENCE 

AND 

SURFACE  CONDITIONS 

!■» 

■DB 

8S# 

1 

80° 

I 

75° 

70° 

I 

FLAT 

3-4 

}flat 

in 

FLAT 

E9 

1 

mm 

em 

FLAT 

3-4 

i 

1  « 

|  0.  804 

0.  110 

|  0.  513 

0.  138 

0.292 

0.  166 

L 

0.  180 

1 

0.  151 1 

0.  120 

0.  120 

KB3I 

jo.  19610. 890 

0.  487 

0.  862 

0.  708 

0.  854 

0.  820 

0. 349 1 

0.  88C 

0.  880 

Ts 

i  12.  41 

12.  10 

12.  72 

13.07 

1 

12.  93 

13.05 

13.04 

13.  06  { 

l 

113.09 

13.  09 

Td 

0.  79 

0.  10 

0.  48 

0.  13 

0.  27 

0.  15 

0.  16 

EDI 

0.  11 

0.  11 

Lz? 

13. 

18  j 

I  Tntd  ! 

0. 

02 

TEMPER. 

EOT  SPOT 

SEA  SURFACE 

AT 

RATIO  | 

;  AGA 

15.  30 

9.  26 

6.  04 

1.  63 

ACTUAL 

16.  89 

13.20 

3.  69 

R  :  REFLECTANCE  OF  THE  SEA  SURFACE 
FLAT:  FLAT  SEA  SURFACE 

8-4  :  SEA  SURFACE  ROUGHENED  3Y  BEAUFORT  4  WIND 

Ts  :•  CALCULATED  APPARENT  SEA  SURFACE  TEMPERATURE 

Td  :  DIFFERENCE  BETWEEN  ACTUAL  AND  CALCULATED  SEA  TEMPERATURES 

Tra  :  AGA  SEA  TEMPERATURE  MEASUREMENT  WITHOUT  ATMOSPHERIC  ATTENUATION 

Trad  :  DIFFERENCE  3ETWEEN  ACTUAL  AND  AGA  SEA  TEMPERATURES 


*  ALL  TEMPERATURES  ARE  IN  °C 
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TABLE  23 

SEA  SURFACE  TEMPERATURE  AND  EFFECTIVE  AT  "E" 


DATS  /  TIMS  : 

WIND  SPEED  /  DIRECTION  : 
AIR  TEMPERATURE  : 

AMBIENT  TEMPERATURE  : 
SEA  SURFACE  TEMPERATURE  : 
CLEAR,  VISIBILITY  10  NM 
ATMOSPHERIC  TRANSMITTANCE 


22  MAY  00:48  -  01:03 
9.5  a/ sec  (3-3)  290° 

9.  7«C 
17 'C 

10.  3“C 

0.  3232 


E  CONDITIONS 

■  ) 

89 4  I 

!  asa 

SO4 

1  7S‘ 

j  70° 

i 

l 

j  FLAT 

3-4  | 

j  FLAT 

3-4  j { FLAT  3-4 

FLAT  3-4 

j  FLAT  3-4 

i  a  ! 

1 

1  0.  304 

o.  noj 

|o. 513 

0.1331  0.  292  !  0.  166  S 

Jo.  ISoio. 151 

0.  120  jo.  120 

|i-a  i 

0.  196 

0.  890 

0.  487 

11  J  i 

0.  362|  |0.  708  j  0.  334  j 

jo.  320 Jo. 349 j 

jo. 380 |0. 380 

Ts 

1  9.  94 1 10.  59  j 

jlO. 28 

10.  67 

JlO. 51 IlC.  64 1 

|1C.  62 1 10.  55 

10.  58 !  10.  68 

Td  j 

|  0.86 

0.  Ill 

{  0.  52  J  0. 13  j  {  0.291  0.  16  j 

|  0.  is!  0.15 

0.  12  J  0.  12 

Tm 

j  10.  76  . 

1 

Tmd 

0.  04 

j  TEMPER. 

HOT  SPOT 

SEA  SURFACE  J 

AT 

RATIO  | 

j  AGA 

8.34 

6. 80  | 

1.  74 

1.  70  1 

J  ACTUAL 

11.  82 

10.  30  1 

1.  02 

r.T.d 


REFLECTANCE  OF  THE  SEA  SURFACE 
FLAT  SEA  SURFACE 

REA  SURFACE  ROUGHENED  BY  3EAUFORT  4  WIND 

CALCULATED  APPARENT  SEA  SURFACE  TEMPERATURE 

DIFFERENCE  3STWEEN  ACTUAL  AND  CALCULATED  SEA  TEMPERATURES 

AGA'  SRA  TEMPERATURE  MEASUREMENT  WITHOUT  ATMOSPHERIC  ATTENUATION 

DIFFERENCE  3STWEEN  ACTUAL  AND  AGA  SEA  TEMPERATURES 


*  ALL  TEMPERATURES  ARE  IN  °C 
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QLLVH  Id 


IX.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

Temperature  and  radiance  distributions  of  an  instrumented  target  measured 
remotely  by  radiometric  imaging  with  the  AGA  Thermovision  780  system  have  been 
compared  with  locally  measured  environmental  parameters..  The  ship-to-background 
temperature  difference,  frequently  used  as  a  measure  of  the  contrast,  was  compared  as 
measured  thermometrically  with  local  temperature  sensors  and  radiometrically  using 
the  Thermovision.  Significant  differences  have  been  observed,  depending  on  the  time  of 
the  day,  windspeed,  and  atmospheric  conditions.  These  differences  have  been  assigned 
to  effects  of  the  low  emissivity  and  high  reflectivity  of  the  sea  surface  at  near  grazing 
angles. 

Reflected  thermal  radiation  from  the  sky  contributed  greatly  to  the  sea  surface 
radiance  sensed  by  the  system.  Comparison  of  mrget-to-background  temperature 
differences  from  radiance  and  from  thermometric  measurements  showed  large 
discrepancies  due  to  the  effect  of  the  reflected  sky  radiation,  characteristic  of  the 
atmospheric  temperature.  This  reflected  radiation  was  strongly  dependent  on  the 
roughness  of  the  sea  surface  under  conditions  of  varying  wind  speed.  The  differences  in 
the  air  and  sea  surface  temperatures  caused  the  effective  (radiometric)  temperature 
differences  between  the  target  and  background  to  deviate  both  positively  and  negatively 
from  the  true  temperature  difference.  With  the  sky  colder  than  the  sea  surface,  as  was 
always  found  in  the  first  series  of  measurements,  the  reflected  sky  radiation  caused  the 
sea  surface  temperature  measured  by  the  Thermovision  to  be  lower  than  the  actual.  As 
a  result,  the  effective  temperature  difference  as  observed  with  the  Thermovision  was 
higher  than  the  actual  difference.  Beth  the  temperature  decrease  and  the  increase  in  the 
effective  temperature  difference  were  related  inversely  to  the  sea  state  conditions. 

An  expression  for  the  sea  surface  radiance  including  the  reflected  sky  radiation 
has  been  developed  and  used  to  model  the  radiometric  effective  ship-to-background 
temperature  difference.  This  model  can  be  evaluated  for  known  atmospheric  conditions 
using  the  LOWTRAN  6  computer  code  and  has  been  found  to  be  consistent  with 
observations  to  date. 
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In  analysis  of  radiometric  images  of  the  ship  taken  at  different  aspect  angles, 
apparent  differences  in  radiance  distribution  were  observed.  This  is  considered  to  result 
from  thermal  radiation  from  hotter  portions  of  the  ship  reflecting  from  cooler  regions, 
increasing  the  apparent  surface  temperature.  Consequently,  temperature  values  from 
only  those  surfaces  more  nearly  perpendicular  to  the  line  of  sight  were  averaged  to 
obtain  representative  temperature  and  radiance  distributions. 

Although  the  quantitative  results  cannot  be  considered  to  be  of  great  accuracy, 
they  show  clearly  the  strong  dependence  of  the  reflected  sky  thermal  radiation  on  the 
•sea  roughness  due  to  the  wind  speed  at  the  times  of  the  observations. 

B.  RECOMMENDATIONS 

Since  we  had  to  calibrate  our  system  empirically  to  achieve  the  accuracy  of  our 
measurements  it  is  obvious  that  we  have  inserted  some  error  factors  into  our  results. 
Therefore,  we  recommend  a  new  series  of  target  and  background  measurements  to  be 
taken  with  the  whole  system  being  accurately  calibrated. 

A  few  months  before  this  work  was  completed  the  Thermovision  was  tested  and 
calibrated  accurately.  The  gain  of  the  preamplifier  in  the  scanner  was  adjusted  and  new 
calibration  curves  were  derived.  However  there  remained  a  problem  in  that  the 
computer  did  not  calculate  the  exact  temperatures  because  the  thermal  values  inserted 
by  the  system  into  the  computer  program  were  off  by  four  thermal  units,  even  though 
a  few'  weeks  before  this  work  was  completed  we  took  a  new  series  of  measurements  of 
the  same  target.  Our  results  and  analysis  are  shown  in  Appendix  A. 

All  our  measurements  were  taken  at  an  incidence  angle  near  the  horizon  where 
the  curves  in  Fig  3.S,  used  to  find  the  reflectance,  are  not  strictly  applicable.  Therefore, 
depending  on  the  operational  use,  we  may  need  this  information  with  greater  accuracy, 
we  recommend: 

For  systems  like  missiles  or  airborne  forward  looking  infrared  devices  a  series  of 
measurements  should  be  taken  for  incidence  angles  between  70°  and  80°  where  the 
curves  in  Fig  3.8  are  applicable,  and  where  more  accurate  results  can  be  obtained. 

For  forward  looking  devices  operating  on  ships,  with  incidence  angles  near  the 
horizon,  long  series  of  measurements,  under  all  possible  kinds  of  environmental 
conditions,  may  be  needed  to  provide  sufficiently  accurate  information. 
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APPENDIX  A 

DATA  ANALYSIS  OF  LATER  MEASUREMENTS 


1.  GENERAL 

A  second  set  of  thermal  imaging  data  were  collected  in  mid  November,  1986.  All 
measurements  were  taken  in  the  spectral  region  8-14  pm  as  before,  using  the  AGA 
Thermovision  supplemented  with  a  microcomputer.  The  ship  target  was  again  the  R.  V 
"POINT  SUR",  scheduled  on  an  operational  oceanography  student  cruise,  and  the 
geographic  location  where  the  thermal  imaging  device  and  the  ship  were  positioned 
during  the  experiment  was  also  the  same. 

These  data  were  collected  in  order  to  obtain  further,  more  accurate, 
measurements,  under  different  environmental  conditions,  after  the  Thermovision  had 
been  calibrated.  This  was  done  to  check  our  previous  results  and  improve  the 
information  on  the  reflected  sky  thermal  radiation  and  the  effective  target-to- 
background  temperature  difference  relation. 

2.  CALIBRATION  TEST 

The  Thermovision  was  accurately  calibrated  after  the  first  measurements  had 
been  taken.  The  preamplifier  gain  of  the  scanning  device  which  was  found  to  be 
incorrectly  set,  was  adjusted  and  new  calibration  constants  were  derived  and  inserted  in 
the  computer  program.  However,  even  if  the  Thermovision  was  operating  accurately  on 
this  occasion,  the  temperatures  calculated  by  the  computer  were  not  accurate.  The 
problem  lay  in  the  thermal  values  transmitted  to  the  computer  from  the  Thermovision. 
The  thermal  values  received  by  the  computer  were  about  3  to  4  thermal  units  lower 
than  measured  directly  on  the  Thermovision.  The  error  was  minimized  by  adjusting 
the  THERMAL  LEVEL  control  to  read  about  4  thermal  units  higher,  using  the 
THERMAL  LEVEL  ADJ  control  on  the  Black,  White  monitor  chassis.  In  this  way  the 
the  computer  received  almost  exactly  the  thermal  values  corresponding  to  the  object 
in  view.  Of  course  this  was  not  the  right  way  to  sohe  the  problem,  since  now  the 
Thermo vi  don  was  measuring  thermal  values  4  thermal  units  higher  than  previously. 
However,  the  system  was  empirically  calibrated  again,  following  exactly  the  same 
procedure  described  in  Chapter  VI. 
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The  calibration  temperature  measurements  for  the  two  thermal  ranges  used  in 
these  series  of  measurements  (  2  and  5  )  am  shown  in  the  following  Tables  25  and  26. 
The  data  were  plotted  as  shown  in  the  following  Figures  A.l  and  A.?  and  the  best  fit 
straight  line  was  found.  As  can  be  seen  from  the  mean  and  variance  of  the  difiVrence 
between  the  exact  and  calculated  temperatures,  the  temperatures  the  computer 
calculated  were  very  close  to  the  eyac:  temperatures,  The  correction  er-or  found  this 
time  was  within  2.9  %  for  the  theimal  range  2  and  5.5  %  for  the  thermal  rang?  5.  Cf 
course  again  we  ha^e  to  keep  in  mind  that  these  calibration  functions  are  valid  on’-y  lor 
temperatures  in  the  region  from  5°C  to  25  cC 

3.  MEASURED  DATA. 

Thermal  images  of  the  target  and  its  background  were  collected  on  the  14^  of 
November  (  09:07  -  09:12  ),  on  the  16^  of  November  (  17:58  -  1S:02  )  and  on  the  18^ 
of  November  (08:34  -  08:36  V  Ail  measurements  were  taken  using  the  same  equipment 
as  before,  in  the  same  spectral  region  and  geographic  location  as  b  described  in 
Chapter  VII.  The  only  difference  was  that  this  time  the  R,  V  "POINT  SUR"dic\  not 
maneuver  mar, the  buoy  to  show  different  aspects  but  passed  by  the  buoy  on  her  wav 
into  and  out  of  the  Monterey  harbor.  The  distance  between  the  lab  station  and  the 
ship  was  about  650  m  on  the  141*1  and  IS^1  when  the  ship  was  going  out  of  Monterey 
Bay  showing  the  port  side  and  about  750  m  on  the  16^  on  her  way  into  the  harbor 
showing  the  starboard  side. 

The  environmental  conditions  experienced  during  these  nuasurements  are  shown 
in  the  data  presentation  tables  for  each  observation  period  (Tables  27,28,29). 
Generally  we  had  low  speed  winds  blowing  from  the  North  West,  overcast  skies  and 
fog  on  the  last  day.  The  air  temperature  was  always  higher  on  this  occasion  (  within 
0.5°C  to  1.9°C  )  than  the  sea  surface  temperature. 

To  analyze  the  data  the  same  procedure  described  in  Chapter  VIII  was  followed. 
The  apparent  sea  surface  temperature  was  calculated  and  the  effective  target-to- 
background  temperature  difference,  sensed  by  the  Thermovision  under  the  local 
environmental  conditions,  was  found.  These  results  are  presented  in  Tables  27,28,29  at 
the  end  of  this  appendix.  Also  a  plot  of  the  ratios  of  the  effective  to  actual  ATs  versus 
wind  speed  is  presented  in  Figure  A. 3  after  the  tables. 


TABLE  24 


THERMAL  RANGE  2 

NEW  CALIBRATION  TEMPERATURE  DATA 


TEMPERATURES  IN  DEGREES  C 


BLACKBODY 

EXACT 

AGA 

COMPUTED 

DIFFERENCE 

5.0 

4.7 

0.3 

6.0 

5.8 

0.2 

7.0 

6.7 

.  0.3 

8.0 

7.S 

0.2 

9.0 

8.7 

0.3 

10.0 

9.7 

0.3 

11.0 

11.1 

-0.1 

12.0 

12.2 

-0.2 

13.0 

12.8 

0.2 

14.0 

13.7 

0.3 

15.0 

14.9 

0.1 

16.0 

15.9 

0.1 

17.0 

17.2 

-0.2 

18.0 

18.2 

-0.2 

19.0 

18.9 

0.1 

20.0 

20.2 

-0.2 

21.0 

20.9 

0  1 

22.0 

22.0 

0.0 

23.0 

23.0 

0.0 

24.0 

24  l 

-0.1 

25.0 

25.2 

-0.2 

MEAN  DIFFERENCE  :  0.02S6 
STANDARD  DEVIATION  :  0.2305 
VARIANCE  :  0.0506 

BEST  FIT  STRAIGHT  LINE  :  Y  -=  (0.9798)  X  -  0.3628 
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TABLE  25 


THERMAL  RANGE  5 

NEW  CALIBRATION  TEMPERATURE  DATA 


TEMPERATURES  IN  DEGREES  C 


BLACKBODY 

EXACT 

AGA 

COMPUTED 

DIFFERENCE 

5.0 

4.6 

0.4 

6.0 

5.5 

0.5 

7.0 

6.6 

0.4 

S.O 

7.6 

0.4 

9.0 

8.4 

0.6 

10.0 

9.6 

0.4 

11.0 

10.5 

0.5 

12.0 

11.4 

0.6 

13.0 

12.5 

0.5 

14.0 

13.6 

0.4 

15.0 

14.6 

0.4 

16.0 

15.7 

0.3 

17.0 

16.7 

0.3 

18.0 

17.7 

0.3 

19.0 

18.6 

0.4 

20.0 

19.6 

0.4 

21.0 

20.7 

0.3 

22.0 

21.7 

0.3 

23.0 

22.6 

0.4 

24.0 

23.7 

0.3 

25.0 

24.6 

0.4 

MEAN'DIFFERENCE 

0.4047 

STANDARD  DEVIATION 

0.0921 

VARIANCE 

0.0081 

BEST  FIT  STRAIGHT  LINE  :  Y  =  (0.9920)  X  -  0.5217 
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AGA  CALIBRATION  CURVE 
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Figure  A. 'l  Thermal  Range  5  new  calibration  temperature  data. 


4.  DISCUSSION  AND  ANALYSIS  OF  THE  RESULTS 

As  was  expected,  reflected  thermal  radiation  from  the  warmer  sky  affected  the  sea 
background  radiance  under  the  various  sea  state  conditions  due  to  the  wind  speeds 
encountered.  Therefore,  the  apparent  and  the  effective  sea  surface  temperatures  were 
found  to  be  higher  than  the  actual.  As  before,  the  differences  between  the  apparent  and 
actual  sea  surface  temperatures  were  higher  at  incident  angles  near  the  horizontal  and 
at  lower  wind  speeds. 

Since  the  effective  sea  surface  temperature  sensed  by  the  Thermovision  was 
higher  than  the  actual,  the  effective  AT  sensed  by  the  detector  was  lower  than  the 
actual  AT. 

The  major  error  factor  occurred  because  we  did  not  have  an  indication  of  the 
actual  temperature  of  the  stack,  the  hottest  spot  on  the  ship,  and  we  had  to  trust  again' 
the  temperature  the  system  calculated  for  this  spot.  On  this  exercise  we  did  not  have 
thermocouples  positioned  at  various  locations  on  the  ship  as  we  had  before,  to 
calibrate  the  measurements.  An  attempt  was  made  to  record  a  measure  of  the  stack 
temperature,  using  a  common  thermometer  fastened  to  the  stack  '  surface. 
Unfortunately  since  we  could  not  place  the  thermometer  on  the  surface  viewed  by  the 
system,  we  did  not  get  accurate  measurements. 

Even  with  all  the  uncertainties  about  the  exact  ATs  and  their  ratios,  the  results 
show  again  that  the  variations  in  emissivity  and  reflectance  of  the  sea  surface  under  the 
local  environmental  conditions  had  affected  the  effective  AT  sensed  by  the 
Thermovision. 

From  the  tabulated  and  plotted  ratios  (Fig  A. 3)  of  the  effective  and  actual  ATs  it 
can  be  seen  that  for  the  lower  wind  speed  of  1.5  m.'s  (  161*1  Nov.),  and  the  higher 
reflectance,  we  find  the  lower  AT  ratio  of  0.939  .  This  ratio  increases  to  0.954  when- 
the  wind  speed  increases  to  2  m.  s  (  IS1*1  Nov.  )  and  the  reflectance  decreases.  Normally 
one  would  expect  to  have  the  higher  ratio  on  the  14^  of  Novembei  when  the  wind 
speed  was  2.9  m/s,  but  unfortunately  the  measured  data  do  not  seem  to  be  consistent 
for  that  date.  The  temperature  value  recorded  for  the  stack  on  this  day  appears 
inconsistent  with  the  reliably  measured  values  of  air  temperature.  However,  the  ratio  is 
close  to  one,  not  too  far  from  what  it  is  expected. 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 

It  is  noticable  again  in  this  series  of  measurements  that  reflected  sky  radiation 
affected  the  sea  surface  radiance  sensed  by  our  system.  Since  the  sky  was  warmer  than 
the  sea  surface  ,  during  the  time  periods  of  observation,  reflected  sky  thermal  radiation 
caused  the  sea  surface  temperatures  sensed  by  the  Thermovision  to  look  higher  than 
the  actual  sea  surface  temperature.  The  higher  effective  sea  surface  temperatures 
caused  the  effective  ATs  sensed  by  the  Thermovision  to  be  lower  than  the  actual  ATs 
under  the  various  wind  speeds  experienced. 

To  improve  the  accuracy  of  the  data  set.  it  is  recommended  that  a  series  of 
measurements  be  taken,  with  the  whole  system  correctly  calibrated  and  with  a 
measurement  of  the  actual  temperature  of  the  side  of  the  stack  viewed  by  the  system. 
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TABLE  26 

SEA  SURFACE  TEMPERATURE  AND  EFFECTIVE  AT  "F 


DATS  /  TIKE 

WIND  SPEED  /  DIRECTION 
AIR  TEMPERATURE 
AMBIENT  TEMPERATURE 
SEA  SURFACE  TEMPERATURE 
TARGET  DISTANCE 


14  NOV  OS:  07  -  OS:  12 

2.9  o/sec  (B-2)  320* 

13. S°C 

15  6C 
12. 2°C 
6  5  Off. 


OVERCAST,  VISIBILITY  4  NM 
ATMOSPHERIC  TRANSMITTANCE  :  0. 8864 


-  -  —  -  -  — l 

ANGLE  OF  INCIDENCE  AND  SURFACE  CONDITIONS 

i 

88° 

!  e5° 

■m 

70' 

|flat 

B-4 

Iflat 

B-4 

Jflat 

B-4 

Jflat 

B-4 

FLAT  B-4 

m 

0.  804 

0.  110 

0.513 

0.292 

0.  166 

Jo. 180 

0.  151 

0.  120  jo.  120 

l-R 

0.  196 Jo.  890 

0.  487 

0.  862 

0.  708 

0.  834 

j  0.  820 

0.  849 

0.  880 jo.  880 

Ts 

m 

iisg 

12.39 

12.  71 

12.  45 

12.48 

12.  42 

12. 36)12.36 

Td 

IB 

am 

DEB 

j-0.  51 

-0.25 

J-0. 28 

-C.22 

j-0. 16 j-0. 16 

Tin 

12.56 

Tmc 

mssm 

11 B 

TEMPER. 

HOT  SPOT 

SEA  SURFACE 

AT 

RATIO 

AGA 

15.  00 

11.  13 

3.  87 

1.  04 

ACTUAL 

IS.  90 

12.  20 

3.  70 

R  :  REFLECTANCE  OF  THE  SEA  SURFACE 
FLAT:  FLAT  SEA  SURFACE 

5-4  :  SEA  SURFACE  ROUGHENED  EY  BEAUFORT  4  WIND 

Ts  :  CALCULATED  APPARENT  SEA  SURFACE  TEMPERATURE 

Td  :  DIFFERENCE  BETWEEN  ACTUAL  AND  CALCULATED  SEA  TEMPERATURES 

Tff  :  AGA  SEA  TEMPERATURE  MEASUREMENT  WITHOUT  ATMOSPHERIC  ATTENUATION 

Tmd  :  DIFFERENCE  3ETWEEN  ACTUAL  AND  AGA  SEA  TEMPERATURES 


*  ALL  TEMPERATURES  ARE  IN  *C 
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TABLE  27 

SEA  SURFACE  TEMPERATURE  AND  EFFECTIVE  AT  "G" 


DATS  /  TIME 
WIND  SPEED  /  DIRECTION 
AIR  TEMPERATURE 
AMSIENT  TEMPERATURE 
SEA  SURFACE  TEMPERATURE  : 

target  Distance  : 


:  16  NOV  17:58  -  18:02 
:  1.5  m/sec  (5-1)  320* 

15. 9*C 
17. S'c 
14°C 
750n 


PARTIALLY  CLOUDY,  VISIBILITY  5  NM 
ATMOSPHERIC  TRANSMITTANCE  :  0. 8842 


i  ANGLE  OF  INCIDENCE  AND  SURFACE  CONDITIONS 


!  85'  | 

mem 

mm 

B 

70' 

(flat 

HE 1 

FLAT 

B-4 

|flat 

B-4  j  j FLAT 

FLAT 

B-4 

j  0.  804 

0.  110 

0.513 

0.  138 

jo. 292 

0.  166 

0.  180 

0.  151 

0.  120 

0.  120 

1-R 

jo. 1S6 

0.  890  j 

0.  487 

0.862 

0.  708 

0.  834 1 

0.  820 

0.  849 

j  0.  880 

0.  580 

Ts 

15.  58 

14.  25 

15.06 

14.31 

14.  63 

14.  74 

14.40 

14.  39  j 

14.27 

14.  27 

Td 

-1.58 

-0.25| 

-1.  06 

-C.  33.  j 

-0.  63  j  -0.  74 

i-0.  40 

-0.39 

-0.27 

-0.27 

Tm  j 

14.95 

Tmd  | 

1— 1 

jl  '  ' 

TEMPER. 

HOT  SPOT 

SEA  SURFACE 

AT 

RATIO 

AGA 

17.38 

13.  22 

4.  16 

0.  939 

ACTUAL 

18.  38 

13.  95 

4.  43 

R  v  REFLECTANCE  GF  THE  SEA  SURFACE 
FLAT:  FLAT  SEA  SURFACE 

B-4  :  SEA  SURFACE  ROUGHENED  BY  BEAUFORT  4  WIND 

Ts  :  CALCULATED  APPARENT  SEA  SURFACE  TEMPERATURE 

Td  :  DIF; -PENCE  BETWEEN  ACTUAL  AND  CALCULATED  SEA  TEMPERATURES 

Tn  :  AGA  SEA  TEMPERATURE  MEASUREMENT  WITHOUT  ATMOSPHERIC  ATTENUATION 

Trod  :  DIFFERENCE  BETWEEN  ACTUAL  AND  AGA  SEA  TEMPERATURES 


*  ALL  TEMPERATURES  ARE  IN  'c 
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TABLE  28 

SEA  SURFACE  TEMPERATURE  AND  EFFECTIVE  AT  "H 


DATS  /  TIME 

WIND  SPEED  /  DIRECTION 
AIR  TEMPERATURE 
AMBIENT  TEMPERATURE 
SEA  SUREACE  TEMPERATURE 
TARGET  DISTANCE 
OVERCAST,  FOGGY,  VISIBILITY  1  NM 
ATMOSPHERIC  TRANSMITTANCE  :  0. 3890 


18  NOV  08:34  -  06:36 

2.0  m/sec  (B-2)  335* 

14.  8*0 

16°  C 

14. 3S°C 

650m 


ANGLE  OF  INCIDENCE  AND  SURFACE  CONDITIONS 


! 

i  89  *• 

85* 

!  80*  I 

75* 

!  70' 

s 

i 

FLAT 

E9 

1 

I  FLAT 

B-4 

Jflat 

B-4  | 

FLAT 

B-4 

FLAT 

ESI 

!  r 

0.  804 

0.  110 

0.  513 

0.  138 

0.292 

0.  366  j 

0.  180 

0.  151 

0.  120 

0.  120 

jl-R  J 

0.  196 

0.  890 

0.487 

0.  862 

0.  708 

0.  834 ! 

0.  820  Jo.  849 

i  0.  880 

0.  880 

Ts 

!  14.  72 

OBI 

mm 

14.  45 

14.  46  j  14.  44 

14.  43  1 14.  43  j 

Td 

j-0.34 

-0.  04 

j-0.22 

-0.06 

1  -0. 12 

-0.071 

j-0.  08 

-0.  06 

-0.05! -0.05 

Tro 

14.  62 

• 

Tmd 

-0.24 

|fl 

TEMPER. 

HOT  SPOT 

SEA  SURFACE 

AT 

RATIO 

AGA 

14.  67 

13.00 

1. 67  J  0. 954 

ACTUAL 

16.  13 

14.  36 

mam 

R  :  REFLECTANCE  OF  THE  SEA  SURFACE 
FLAT:  FLAT  SEA  SURFACE 

B-4  :  SEA  SURFACE  ROUGHENED  BY  BEAUFORT  4  WIND 

Ts  :  CALCULATED  APPARENT  SEA  SURFACE  TEMPERATURE 

Tc  :  DIFFERENCE  BETWEEN  ACTUAL  AND  CALCULATED  SEA  TEMPERATURES 

Tm  :  AGA  SEA  TEMPERATURE  MEASUREMENT  WITHOUT  ATMOSPHERIC  ATTENUATION 

Trod  :  DIFFERENCE  BETWEEN  ACTUAL  AND  AGA  SEA  TEMPERATURES 


*  ALL  TEMPERATURES  ARE  IN  cC 
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WIND  SPEED  TO  DT  RATIO  RELATION 


oiivh  xa 
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Figure  A. 2  AT  ratios  plotted  versus  wind  speed. 


APPENDIX  B 

IN-BAND  FLUX  CALCULATION  USING  TI-59  PROGRAM 


1.  GENERAL 

The  program  for  the  TI-59  calculator  used  to  calculate  the  in-band  flux  makes 
use  of  the  Simpson's  approximation  program  included  in  the  master  library  module.  It 
was  originally  written  for  the  HP-67  calculator  by  Robert  Pitlak  (  Appolo  Lasers,  6357 
Arizona  Circle,  Los  Angeles  CA.,  90045  ).  and  has  been  translated  for  the  TI-59. 

To  use  this  program  the  following  values  have  to  be  inserted: 

a.  The  initial  and  final  wavelengths  of  the  spectral  band. 

b.  The  temperature  of  the  object  ( in  °K  or  °F  or  °C  ). 

c.  The  even  number  of  intervals  to  be  used  bv  the  Simpson's  approximation  of  the 
integral. 

The  calculated  flux  is  given  in  Watts/cm  or  photons, 'cm'sec 

2.  FORMULAE  AND  CONSTANTS 

As  we  discussed  in  Chapter  VI,  Section  E,  Planck's  blackbody  law  can  be  written 
in  the  form 

X2 

Q  =  J  CQ  q  dX  photons/cm^sec 

x2 

or  M  =  jCj(q.'X)dX  Watts.'cm2 


Where  q  =  -  x  =  - 

X4(ex-1)  XT 

Cq  =  2kc  =  1.883651556  X  102  ( fim2/cm2sec  ) 

Cj  =  2tthc2  =  37418.42875  (  Watts  ftm4.'cm2  ) 

C2  =  ch/k  =  14388.32334  ( *tm  °K  ) 

The  integrations  are  done  numerically  using  Simpson's  approximation. 
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3.  PROCEDURE 

How  to  use  the  program 

a.  Load  the  magnetic  card 

b.  Press  RST 

c.  Enter  Xj  initial  wavelength  in  fim.  press  R/S 

d.  Enter  final  wavelength  in  Urn,  press  R/S 

e.  Enter  temperature  (T).  press 

(1)  A  if  in  0  K 

(2)  B  if  in  0  F 

(3)  C  if  in  °C 

f.  Enter  even  number  of  integration  intervals 

g.  Press 

(1)  D  to  get  the  answer  in  Watts, 'cm^ 

(2)  E  to  get  the  answer  in  photons/ cm~sec 

4.  PROGRAM  LISTING 

In  the  following  Table  24  we  give  the  listing  of  the  program  used  in  the  TI-59 
calculator. 
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TABLE  24 


TI-59 

PROGRAM  LISTING 

114 

02 

2 

174 

'A -A 

l  h:i 

055 

07 

7 

i 

115 

03 

■V 

175 

370 

056 

03 

3 

1 16 

03 

3 

176 

06 

06 

057 

35 

a 

117 

04 

4 

i  *• 

*  i  « 

36 

PCM 

053 

42 

370 

US 

55 

* 

173 

03 

0? 

000 

42 

370 

053 

00 

00 

113 

43 

RCL 

17? 

14 

D 

001 

01 

01 

060 

I1 

R.'S 

120 

07 

07 

130 

22 

IMV 

002 

3° 

PR7 

061 

76 

LSL 

121 

55 

~ 

131 

C*  A 
xi  i- 

EE 

003 

?i 

R.'3 

062 

13 

C 

122 

43 

RCL 

132 

PR7 

004 

42 

370 

063 

42 

370 

123 

00 

00 

133 

33 

AD1-' 

005 

02 

02 

064 

00 

00 

124 

54 

> 

134 

31 

R.'C 

006 

33 

PR.  7 

065 

01 

l 

125 

54 

> 

135 

7  6 

LSL 

oo  r 

31 

R.'3 

066 

05 

5 

126 

75 

- 

136 

15 

£ 

003 

76 

LSL 

067 

00 

0 

127 

01 

i 

137 

33 

ftDV 

003 

11 

fi 

063 

00 

0 

123 

54 

) 

133 

42 

370 

010 

42 

370 

063 

00 

0 

123 

.  54 

) 

IS? 

05 

05 

Oil 

00 

00 

070 

00 

0 

130 

1  ,u 
l  •  n 

130 

•50 

PR7 

012 

02 

071 

63 

OP 

131 

65 

X 

131 

33 

ADV 

013 

06 

6 

072 

04 

04 

132 

43 

RCL 

132 

53 

< 

014 

00 

0 

073 

43 

RCL 

133 

03 

03 

133 

43 

RCL 

015 

nn 

0 

074 

00 

00 

134 

54 

> 

134 

02 

02 

016. 

00 

0 

075 

6? 

OP 

135 

32 

R7M 

135 

75 

- 

017 

00 

0 

076 

06 

06 

136 

7  6 

LSL 

136 

43 

RCL 

013 

63 

OP 

077 

63 

OP 

137 

14 

D 

137 

01 

01 

013 

04 

04 

073 

00 

00 

133 

33 

hDV 

133 

54 

> 

020 

43 

RCL 

073 

43 

RCL 

133 

42 

370 

TOO 

i  «  < 

55 

021 

00 

00 

030 

00 

00 

140 

05 

05 

200 

43 

RCL 

022 

63 

OP 

031 

35 

* 

141 

33 

PR7 

■201 

05 

05 

023 

06 

06 

032 

02 

"A 

142 

ftft 

fiSV 

202 

Q5 

3 

024 

63 

OP 

033 

07 

“T 

143 

53 

( 

203 

42 

370 

025 

00 

00 

Oft  4 

03 

3 

144 

43 

RCL 

204 

03 

03 

026 

31 

R'"3 

035 

35 

a 

145 

02 

02 

205 

01 

1 

027 

76 

LSL 

036 

42 

370 

146 

75 

- 

206 

33 

• 

023 

12 

S 

nsr 

00 

00 

147 

43 

RCL 

207 

03 

3 

023 

42 

3’0 

033 

31 

R/3 

143 

01 

01  • 

203 

03 

3 

030 

00 

00 

033 

7  6 

LSL 

14? 

54 

> 

20? 

03 

3 

031 

02 

A 

030 

16 

fi* 

150 

eir 

- 

210 

04 

6 

022 

01 

I 

031 

42 

370 

151 

43 

RCL 

211 

05 

5 

023 

00 

0 

032 

07 

07 

152 

05 

05 

212 

01 

i 

034 

00 

o 

033 

53 

< 

153 

35 

3 

213 

06 

6 

035 

00 

0 

03  4 

53 

< 

154 

42 

370 

214 

52 

EE 

036 

00 

0 

035 

43 

RCL 

155 

03 

03 

215 

02 

037 

s? 

OP 

036 

07 

07 

156 

03 

216 

03 

3 

033 

1.14 

1.14 

037 

45 

Y“ 

157 

33 

• 

217 

42 

370 

033 

RCL 

033 

43 

RCL 

153 

07 

% 

213 

03 

03 

040 

00 

00 

033 

03 

03 

15? 

04 

4 

21? 

04 

4 

041 

63 

OP 

.  100 

65 

X 

160 

01 

1 

220 

42 

370 

042 

06 

06 

101 

53 

< 

161 

05 

5 

221 

0? 

0? 

043 

6? 

OF 

102 

53 

< 

162 

00 

0 

01 

1 

044 

00 

00 

103 

43 

RCL 

163 

07 

*• 

V 

A  A 

IMV 

0'45 

75 

- 

iru 

06 

06 

164 

04 

4 

224 

23 

LUX 

046 

03 

3 

105 

45 

yx 

1 65 

tr  a 
Vi. 

EE 

225 

42 

370 

047 

02 

*» 

106 

53 

< 

1 66 

04 

4 

226 

06 

06 

043 

35 

a 

107 

01 

1 

167 

42 

370 

227 

36 

PCM 

043 

65 

X 

103 

04 

4 

163 

03 

03 

223 

0? 

09 

050 

05 

5 

103 

03 

3 

163 

05 

5 

AAO 

14 

D 

051 

55 

- 

110 

03 

3 

170 

42 

370 

230 

'A'A 

IMV 

052 

03 

7 

111 

03 

3 

171 

03 

03 

231 

52 

££ 

053 

35 

+ 

112 

33 

172 

01 

1 

AAA 

3? 

PRT 

054 

02 

O 

«• 

113 

03 

3 

173 

IMV 

223 

33 

ADV 

234 

91 

R-'S 
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